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I. INTRODUCTION 


This paper presents the development of a computer simulation of a specific 
example of a Jackson network. 

In a typical computer simulation, a series. of mathematical and logical 
relationships are used to create a model of some real world system. The model is then 
exercised over some period of time, and data are collected to estimate the desired 
characteristics of the real word system. When the simulation run is completed, the 
desired statistics are printed out. 

What makes this simulation different, 1s that it also provides the user with a 
visual picture of what 15 occuring in the queue network as the simulation progresses. It 
accomplishes this by actually displaying on the computer screen the arriving and 
departing customers and showing how they move through the network. This portrayal 
is possible because the simulation gives the observer the impression of real-time 
passage. 

Because Jackson networks, and their associated statistical quantities of interest 
can be solved in closed mathematical form, there 1s little gained in this respect from 
simulation. However, by allowing a student to actually observe the inner workings of 
the queue network as it progresses in time, insights into the process can be obtained 
that are not possible from conventional instructional tools alone. 

This paper is the third known effort at the Naval Postgraduate School to develop 
visual simulations of various stochastic processes utilizing the personal computer. 
Prior undertakings are a simulation of various models that are based on an underlying 
Poisson process [Ref. 1], and a simulation of a Machine-Repairman Model [Ref. 2]. 
These two papers provided many ideas for the writing of the Jackson network 
simulation and served as an outline for this paper. 

The remainder of this paper describes the Jackson network model and its 
implementation on the IBM Personal Computer (IBM PC). Chapter IV is a user's 
guide which contains instructions for actually running the simulation. A program 


listing 1s contained in Appendix A. 


II. THE JACKSON NETWORK MODEL 


A. GENERAL DISCUSSION 


In a network of queues, a customer! 


may receive service from one or more 
service centers, possibly revisiting some or all of the service centers many times before 
finally exiting the network. With no additional assumptions made, this type of queue 
network can be extremely difficult, if not impossible to analyze in closed mathematical 
form. In these cases, computer simulation may be the only alternative [Ref. 3: p. 225]. 

One type of network that does permit analysis without the need for simulation is 
the Jackson network. A Jackson network 1s a collection of queues that possess the 
following characteristics [Ref. 4: p. 456]: 


e The network consists of N service centers numbered 1,2,...,|N. Service center 1 
contains C; identical servers. 


e Customers from outside Nu Du. arrive at service center 1 according to a 
Poisson process with rate a I» dw. 


e After receiving service at center i, a customer leaves the network with probability 
P.o > 0 or goes instantaneously to service center J with probability Pij NT -— 
These probabilities are independent of the history of the network. 

e Customers ΓΙΝΕ at service center i (from inside or outside the network) are 
served according toa First, Come First Serve (FIFO) queue discipline. The 
Service times are independent, identically distributed exponential random 
variables with mean 1/μ., = Oe. 

e There is unlimited waiting space at every service center. 

e The effective arrival rate at each service center taking into account arrivals from 
both outside and inside the network must be less than the corresponding service 
rates. 

Two further characteristics of a Jackson network result from the requirement that 
customers from outside the network arrive at the service centers according to a Poisson 
process. This implies that the arrival process must have independent and stationary 
increments. An arrival process is said to possess independent increments if, for 
example, the number of arrivals that occur between times 5 and 10 are independent of 
the number of arrivals between times 35 and 40. An arrival process has stationary 
increments if the distribution of the number of events which occur in any interval of 


time depends only on the length of the time interval. This is the same as saying that 


“Although this paper will refer to customers, implying people, the flow through a 
queue network could just as easily be an object or a job. 
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there is no time of day effect on the arrival rate of customers. One further assumption 
of a Jackson network is that the servers at the different service centers work 


independently and do not assist one another, even if one 1s idle. 


B. THE MODELED JACKSON NETWORK 

The simulation described in this paper is a specific example of a Jackson network 
τω --- ] a k, = Ou, να 0, P... πα. 
GNE ο O s Poo πη, 5 
this network is provided in Figure 2.1. In words, a customer arriving from outside the 


where N = 2, c 


A pictorial representation of 


network must first go to service center | where he enters a FIFO queue. Upon 
completion of service at service center l, the customer can either exit the system with 
probability P,, or proceed to service center 2 with probability 1 - P,,, where he again 
enters a FIFO queue. Upon completion of service at service center 2, the customer 
can either exit the network with probability P,, or loop back to service center 1 with 
probability 1 - P, again entering the FIFO queue. In this configuration a customer 
can loop from service center | to service center 2 and back to service centenal 


numerous times before exiting the network. Each service center has only one server. 





Pio =x it Probabae rey 
A = Arrival Rate 
P 

One Server One Server 20 
i, = Service u;= Service 

Rate ate 

Service Center l Service Center 2 
Ρο1 = Feedback Probabillty 


Figure 2.1 Jackson Network With Two Service Centers. 
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[t is arguable whether any real world system would ever meet all of the 
requirements of a Jackson network. But then that is always the predicament faced by 
the modeler; making the model represent the real world without making it 
unreasonably complex. There are numerous real world systems that might be 
reasonably modeled by this simple Jackson network. For example, the network might 
describe the flow of a piece of equipment through a repair facility. Upon entering 
Service center 1, the repairman evaluates the piece of equipment and attempts to repair 
it. If it can not be repaired, it leaves the network. If it is repairable, it is repaired and 
proceeds to the Quality Assurance Shop (service center 2). At Quality Assurance, the 
repaired item undergoes testing and is either classified as fit or unfit. If it is determined 
fit, it departs the network and returns to the field. If it is determined unfit, it is 
returned to service center | for additional repair. 

Another possible example would be a computer system. Jobs arrive at the 
central processing unit (service center 1) of a computer. When a job finishes CPU 
service it either departs the computer system or requires Input/Output service (service 
center 2). When a job completes IO service it either leaves the computer system or 
returns to the CPU for additional service. 

For this Jackson network model, the state of the system, N(t) = n, 1s defined as 
the total number of customers in the system both waiting and being served at time t. 
This state variable is further broken down into the total number of customers waiting 


and being served by service center 1, N,(t) = n,, and the total number of customers 


1) 
waiting and being served by service center 2, N,(t) — n, [t follows that 
N(t) = N,(t) + N,(t). Because arrivals and departures are discrete events, the value 
of N(t) will always be zero or a positive integer. However, because there is unlimited 
waiting space in front of the service centers, there is no upper limit on the value that 


N(t) can assume. That is, n can be any integer value between 0 and positive infinity. 


C. ANALYSIS OF JACKSON NETWORKS 

It is not the intent of this section to develop in detail the theoretical analysis of 
Jackson networks. However, a general discussion will be included for completeness. 
The interested reader should refer to the reference list and the bibliography for books 
on queueing theory, paying particular attention to the chapters on networks of queues. 
Of particular interest might be the paper written by J.R. Jackson for whom this type of 


network is named. 
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By imposing the restrictions described in Section I-A on a network, Jackson 
showed that in the steady state each service center (say center 1) in the network behaved 
as if it were an independent M/M/c queue with service rate u, and Poisson input rate 
R;. The input rate, R;, can be considered the effective arrival rate at service center 1. 
The effective arrival rate to a service center 1s the sum of the Poisson arrivals from 
outside the system (exogenous customers), plus the arrivals from any other service 
center in the network. The effective arrival rates can be found by solving the following 


set of simultaneous equations: 


In order to prevent an unstable or explosive queue from occurring at any of the 
service centers, It is required that R, < c;u,. In words, the effective arrival rate of 
customers at service center 1 must be less than the cumulative service rate of the servers 
at service center 1. 

The implications of Jackson’s findings are that a complex network can be 
decomposed into a number of simpler M/M/c sub-systems. These sub-systems can 
then be analyzed using Little’s formula and the standard results for the M/M/c queue. 

An interesting aspect of this result is that it can be proved that allowing feedback 
in the network destroys the Poisson flow of arrivals [Ref. 4: p. 460]. However 
Jackson's results show that even for networks that allow feedback, the service centers 
behave as if they were fed totally by Poisson arrivals, when in fact they are not. 
Statisticians have been unable to provide an explanation as to why this is true [Ref. 5: 
p. 330]. 


D. MODEL ESTIMATIONS 

The purpose of any simulation is to gain some insight into the real world system 
being modeled. This insight is usually gained by estimating certain quantities 
pertaining to the systems performance. As a part of the graphical output, the 
simulation provides a table of estimates for a number of these quantities. Because 
Jackson networks can be analyzed in closed form, the graphical output will also include 
the theoretical values for these same quantities. This will provide the user with a 


measure of how well the computer simulation is modeling the real world system. 
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The following equilibrium performance quantities will be both estimated and 


determined theoretically by the simulation: 


B5 = the effective arrival rate of customers at service center 1(1 —^ 1,2). 


B the expected number of times a customer will pass through service 
center 1 before exiting the network (1 = 1,2). 


L; = the expected number of customers both waiting and being 
served at service center 1 (1— 1, 


L = the expected number of customers in the network, both waiting 
and being served. 


W: — the expected amount of time a, customer spends waiting and 
pemg served at service center 1(1 —- |, 

W = the total expected amount of time a customer spends in the network, both 
waiting and being served. 

Z = the proportion of time service center one is busy and Service 


center two ıs idle. 


The theoretical values will be represented by the appropriate capital letter. The 
estimated values will be represented by the appropriate capital letter with a hat over it. 

The following quantities are parameters of the model which can be selected by 
the user prior to running the simulation and are used in determining the theoretical 
quantities just listed. The parameter names are consistent with those used in the 
simulation model. 


A = the Poisson arrival rate of exogenous customers arriving at service center | 
(previously referred to as 


15 
li; = the exponential service rate for service center (i= 1,2). 
; = the probability of exiting the network after completing service at center 1, 
1=1,2 (previously referred to as P;,). 
| - P, = the probability of proceeding to service center 2 after completing 
service at service center | (previously referred to as P,,). 
| - P, = the probability of proceeding to service center | after completing 


service at service center 2 (previously referred to as P,). 


For use in estimating some of the performance quantities, À and p; will also be 


estimated by the simulation model. 


The values of R, and R,, the theoretical effective arrival rates at service center | 


and service center 2, are found by solving equation 2.1, namely: 


| 


R 


πα... (eqn 2.2) 


R 


,7 0 * ΕΙ(Ι-ΡΙ), (eqn 2.3) 


which results in: 


ve 
| 


1 AICP UP (eqn 2.4) 


"3 
| 


EE (eqn 2.5) 


The estimated value of R, is computed by recording the total number of 
customer arrivals at service center l, both exogenous customers and feedback 
customers from service center 2, and dividing by the total elapsed time that the 
simulation has run. The estimated value of R, 1s calculated in a similar manner except 
that there is only one source of arrivals to service center 2, namely those customers 
departing service center | who do not exit the network. 

The estimated value of À, the Poisson arrival rate of exogenous customers 
arriving at service center I is computed by recording the number of exogenous 
customers arriving at service center | and dividing by the total elapsed time that the 
simulation has run. 

The estimated value of p, the exponential service rate for service center l is 
computed by recording the total number of customers completing service at service 
center l, and dividing by the total amount of time service center 1 was busy ( NOT the 
total elapsed time of the simulation). The estimated value of p, is computed in a 
similar fashion. 

The value of L, and L,, the theoretical expected number of customers either 
waiting or being served at service center 1 and 2, are found by utilizing Little’s law and 


the results for a M/M/1 queue: 


a n (eqn 2.6) 


E o OB). (eqn 2.7) 


It follows that L, the expected number of total customers in the network, both 


waiting and being served is simply the sum of L, and L,: 
L D UL. (eqn 2-9) 


. A A. A | j A A ^ ^. 
The estimates L, ,L, and L are found by substituting R,, R,, pi, and ji, into the 


above equations to obtain: 


A A An 

L. M R,/(8,-R.) (eqn 2.9) 
A A A A 

EUR K.-R,) (eqn 2.10) 
^ ^ ^ 

Em. tL, (eqn 2.11) 


The values of W, and W,, the theoretical expected time a customer spends 
waiting and being served at service center | and 2, are also found by utilizing Little’s 


law and the results fora M/M/1 queue: 


Be —- L/R 


m ἢ (eqn 2.12) 


(eqn 2.13) 


The value of W, the total expected time a customer spends in the network, both 


waiting and being served is also computed using Little’s formula: 


W= LA. (eqn 2.14) 
I A A ^ ο NSN AN Δ 
The estimates W,, D and W are found by substituting L, Ex E Ri» Ro and À 
into the above equations to obtain: 
A _ A A 
Do - l; R, (eqn 2.15) 


ο 


A A A 

Wo R (eqn 2.16) 
A AA 

W = LA. (eqn 2.17) 


The value of Z, the theoretical long run proportion of time that service center | 
is busy and service center 2 is idle, is found by using a result commonly referred to as 
Jackson’s Theorem. This theorem, when simplified for the situation where each service 
center has exactly one server, states that for a Jackson network with R; <u; and with k 


service centers the limiting probability is given by: 


Prob(n; at service center 1, 1=1,2,...k) = m(n,,n,,...n,) = (eqn 2.18) 


π(πρπίηρ)..πίηρ) - Π(Βγμρ᾽χ!-Βγμῳ. 


Where n; is the number of customers waiting and being served at service center 1, and 
™(n;) 1s the long run probability of finding n; customers at service center i. The 
theorem says that the limiting probability is equal to the product of the marginal 
probabilities for the individual M/M/1 queues comprising the network. 

To determine the value of Z, it 1s necessary to determine the probability that 
there are no customers at service center 2 and one or more at service center 1. This 


probability is determined as follows: 


Zo Aes) O(n) (eqn 2.19) 
= (L(R 0 O (R I O RO Penis (eqn 2.20) 
= M m. (eqn 2.21) 

a eee ae | 
The estimate Z 15 found by substituting Κι, Κο, μι, and H, into equation 2.21 to 
obtain: 
νο ο. A 
A μοι P qn.) (eqn 2.22) 
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The value of S,, the theoretical expected number of times a customer will pass 
through service center | before exiting the network is found by recognizing that the 
underlying distribution is geometric. For service center 1, it can be seen that the 
following relationship holds, where N, equals the number of times a customer visits 


service center | before exiting: 


ee — IP, * (I-P)P,J (egn 2.23) 
P(N, =2) = {P, + (1-P))P,}{(1-P,)(1-P,)}° (eqn 2.24) 
P(N,=3) = (P, + (1-P,)P,)((1-P,X1-P,)? (eqn 2.25) 
EO =n) = (P + (1-P,)P,)((1-P,NI-P,))" , n= 1,2,3.... (eqn 2.26) 


Expandinzg the right hand side of equation 2.26, 
E - 1 n-1 
ο ο (Po +P -P oP.) (IP -P.+P; P.) (eqn 2.27) 


Letüng 0 = 1-P, -P, + P,P,, it follows that the quantity raised to the first power in 


equation 2.27 is equal to 1 - 0. Equation 2.26 can be rewritten in the form: 
P(N=n) = (1-0(0)??, n = 1,2,3.... (eqn 2.28) 
Which is recognized as a geometric distribution with an expected value equal to: 


"= Do) mE t Po EE): (eqn 2.29) 
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A similar type of analysis can be done to find the theoretical value of S,. Letting 
N, equal the number of times a customer visits service center 2 before exiting and 


8 = (1-P,)(1-P,) as before, it can be seen that the following relationships hold: 


P(N, =0) = P, (eqn 2.30) 
P(N, =1) = OP, + (1-P,)P, (eqn 2.31) 
P(N 2) 205p T O ERR (eqn 2.32) 


n — nn H-1 
PON, = 0) ss Ε.Ε ἪΝ (eqn 236) 
Taking the expected value: 
E(N,) = 0P, ἜΘ; ΟΡ το (eqn 2.34) 


+ 1(1-P,)P, + 20/(1-P,)P, + 30°(1-P,)P, + ...+ nO®!(1-P,)P 
UP» pP» pP» pP» 


ο. + 201530 m n, (eqn 2.35) 
* (POPE 20--6:30- s n, 


Noting that the geometric sequences converge in the limit to (1/1-0)^, equation 2.35 


can be rewritten as: 
E(N,) = P,6(1/1-6)* + (1-P,)(P,)(1/1-6)?. (eqn 2.36) 
The final result after manipulation of equation 2.36 is: 


S, = E(N,) = (1-P,)/(P,+P,-P,P,) = (1-P.)S,. (eqn 2.37) 


The estimates, Ei and $a are found by using the definition of expected value for 
a discrete random variable. Let P:(n) be a shorter notation for the probability that a 
customer is served exactly n times at service center 1 before he departs the network. 
The estimate Pin), Is determined by recording how many customers were served 
exactly n times by service center 1 before departing the network and then dividing this 
number by the total number of customers who have departed the network. The 


estimates of S, and S, are then determined by: 


^ OO A 
S, = 2nP, (n) (eqn 2.38) 
^ OO ^ 
S, = 2nP,(n) (eqn 2.39) 


It should be noted that the simulation program, JACKQUE, is Incapable of 
explicitly tracking an individual customer for more than nine visits to either service 
center. All customers who exit the network having visited either service center, ten or 
more times, are counted as if they had visited that center exactly ten times. This 
results in $, and S, being under-estimates of the theoretical values. Section III will 
discuss the implications of this model limitation. 

Whenever a quantity is estimated through simulation, the estimate will be a point 
estimate of the true value. This point estimate, being a random variable, will have 
some variation associated with it. As the simulation proceeds in time, the point 
estimate will fluctuate around the true value. At any particular point in time the point 
estimate may or may not be close to the true value. To overcome this, a simulation 
model 1s typically run a large number of times using some stopping rule. The ending 
estimates for each run are then averaged to get more accurate estimates. Some sort of 
variance reduction technique like Jack-Knifing or Boot-Strapping is also typically 
employed. 

In JACKQUE however, the simulation model is exercised only once, then 
terminated. For this reason, at any particular point in time, the estimates might not 
coincide closely with the theoretical values. This is particularly true of the estimates 
for L, W, and Z which are algebraic combinations of other point estimates. For 
example, the quantity W, is estimated using the quotient T/R.. If at the point in time 
the estimate 1s made, jm is larger and R is smaller than their true values, the resulting 


estimate of Z will be very much larger than it's true value. 
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Since the primary purpose of JACKQUE is not estimation, this is not a serious 


shortcoming of the model. 
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III. THE SIMULATION 


KWA GENERAL 

The Jackson network is modeled as a discrete-event simulation model. The 
simulation is both dynamic and stochastic in nature. A discrete system (vice 
continuous) is one in which the state variables can only change at a countable or finite 
number of points in time. These points in time are the ones at which an event occurs. 
For the Jackson network, the state of the system N(t), can only change as a result of 
an exogenous arrival, or a service completion by either service center. 

A dynamic simulation (vice static) is one that changes or evolves in time. The 
Jackson network, which is driven by inter-arrival times and service times, 1s obviously 
of this nature. 

A stochastic simulation model (vice deterministic) is one that contains one or 
more random variables. In other words, some characteristic of the system is random in 
nature. For the Jackson network, both the exogenous customer arrival rate, and the 
Service rates are random variables generated from an underlying exponential 
probability distribution. 

A dynamic simulation requires some mechanism for Keeping track of the passage 
of elapsed time, allowing the program to know when, and in what order events will 
occur. In the case of the Jackson network, we need to keep track of when arrivals and 
service completions occur. The variable in the simulation model that performs this 
task is called the simulation clock. The unit of time for the clock is usually not 
explicitly stated, but needs to be consistent with the units of the input parameters. 

There are two basic approaches for advancing the simulation clock. It can be 
accomplished using a next-event time advance, or a fixed increment time advance. In 
both cases the clock is initialized to zero and the times of all future events are 
determined. 

With a fixed-increment clock, the clock 1s advanced in increments of At. After 
each At increment, a check is made to determine if any events should have occured. If 
One or more events should have occured, these events are considered to have occurred 
at the end of the interval, and the system state, N(t), is updated accordingly. Two 


obvious problems with a fixed-increment clock are that event times appear to be larger 
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than they should be due to the rounding up effect, and if more than one event occurs 
in At it might not be possible to decide in which order the events occurred. 

With a next-event clock, the successive jumps of the clock are unequal in size, 
and correspond to the amount of time till the next event occurs. After the simulation 
clock is advanced to the first of these future events, the state of the system N(t) is 
updated to account for the occurrence of the event, and times of future events are 
generated as necessary. The next-event clock is used by most simulation languages 
because it does not introduce bias into the estimates, and by skipping over periods of 
inactivity the program runs faster. 

For purposes of the graphical display, JACKQUE uses a fixed-increment clock. 
The clock is increased in increments of At = 1. This gives the viewer of the graphics 
the illusion of real time passage. This type of clock is also necessary so that the system 
state-versus-time graphs can be updated with fixed length horizontal line increments. 
Although the clock 1s incremented in fixed increments, and all events are plotted as if 
they occurred at the end of At, the program does preserve the correct chronological 
order of the events when plotting. 

For purposes of the estimated performance measures, the program keeps track of 
continuous time, and does not round event occurrences up to At. In this sense, the 


program behaves as if it is using a next-event clock. 


B. THE COMPUTER 

The computer used for this simulation is the IBM Personal Computer (IBM PC), 
equipped with a color monitor. In order to run the simulation the IBM PC must be 
configured with at least 128 kilo-bytes of random access memory (RAM), one 5 1/4” 
floppy disk drive, and a Color/Graphics adapter. The color monitor and 
Color/Graphics adapter are essential since the program uses color graphics. 

The simulation program, JACKQUE was written using version 2.1 of the 
PC-DOS disk operating system. Any version of PC-DOS earlier than 2.1 will also 
work. Because JACKQUE directly interacts with the disk operating system, there can 


be no guarantee that the simulation will run on IBM Compatible machines. 


C. THE PROGRAMMING LANGUAGE 
JACKQUE is written in IBM Advanced Basic (BASICA). This language is 
provided with PC-DOS and possesses all of the graphics capabilities necessary for the 


simulation. This version of Basic is an interpreter language and runs slower than a 
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compiled language.2 However since the illusion of elapsed time is desired in this 
particular application, slower program run time 1s not a disadvantage. 
JACKQUE also uses one sub-routine, called RNGEN, which 1s written in 8088 


assembly language. 


D. PROGRAM STRUCTURE 

The program listing for JACKQUE is contained in Appendix A. The program, as 
it appears in the Appendix, has been extensively edited with comments to assist the 
interested reader in understanding the program. The program contained on the 
distribution diskette does not include the comments, since the inclusion of the 
comments exceeds the memory available within the Basic stack. Additionally, extensive 
use of subroutines were used in an effort to facilitate following the program flow. The 
program consists of five major modules. These modules are labeled Main Program 
Module (MPM), Clock Module (CM), Arrival Module (AM), Service Center 1 Module 
(SCIM), and Service Center 2 Module (SC2M). Any subroutine called exclusively by 
one of the five modules is physically located within that module. After the fifth 
module, all subroutines that are called from more than one module are listed, and 
labeled accordingly. All subroutines have been labeled with appropriate titles. After 
each subroutine title there is a slash (/), followed by the abbreviation of the module or 


modules from which it is called. 


E. DESCRIPTION OF GRAPHICS DISPLAY 

There are two configurations of graphics available to the user. When the 
simulation 1s first started, the graphics display screen will appear as seen in Figure 3.1. 
The upper half of the screen contains a pictorial representation of the status of the 
Jackson network for a particular point in time. A circle in a box represents a person 
being served at that service center. A circle in front of a box represents a person 
waiting at that service center. There may be more customers waiting than can be 
physically displayed on the screen. The number contained in the circle indicates the 
number of times that the particular customer has looped back to service center | from 
service center 2. The absence of a number indicates that the customer is just entering 


the network and has not yet experienced any feedback. 


Ἢ . . . o . 
^There are compiled versions of Basic in existence 
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Figure 3.1 Screen Display with Statistics Table. 


In Figure 3.1, there are twelve customers waiting at service center | as indicated 
by Q, = 12, in addition to the one customer currently receiving service. At service 
center 2 there are four customers waiting for service in addition to the one customer 
receiving service. Of the customers visible, four have experienced feedback twice while 
one has experienced feedback once. 

The bottom left corner of the screen contains the state-versus-time graph for 
service center |. The system state values for the service center are displayed on the 
ordinate of the graph, and 32 increments representing time on the abscissa. The 
display shows that at time = 3818, the state of service center | jumped from 
State 8 to 7. 

The table in the lower right corner contains the estimated and theoretical 
quantities discussed in Section Il-D. The variable names appearing in this table are 
consistent with those defined in Chapter II. 

Once the simulation is running, the user can change the configuration of the 
graphics screen to that shown in Figure 3.2. In this case, the state-versus-time graph 
for service center 2 is displayed in the lower right corner of the screen. If the 
state-versus-time graphs for both service centers were superimposed on each other, the 


resulting graph would be the overall system state-versus-time graph. 
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Figure 3.2 Screen Display with out Statistics Table. 


A more complete descrption of the graphic displays, and user commands 


affecting them, are contained in Chapter IV, the user's manual. 


F. PROGRAM DESCRIPTION 
The main program module contains numerous subroutines that accomplish most 
of the routine housekeeping chores. These include: 


e Checking to see if BASICA has been loaded from DOS, and if the machine is 
configured with a color/graphics adapter. 


e Printing the title screen. 

e Loading user defined functions. 

e Loading the random number generator sub-routine. 

e Initializing and dimensioning variables. 

e Printing the main program menu. 

e Creating symbols that will be used latter in the graphic display. 

e Generating the graphics display. 

e Computing the theoretical values of the tabulated quantities. 

The Program Menu presents the user with a number of options: (1) view the 

program instructions, (2) select a new seed for the random number generator, (3) start 


the simulation using the default parameters, (4) change the model parameters, and (5) 
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cnd the program. The on-screcn instructions provide a condensed version of the user 
instructions found in Chapter IV, the user's manual. The main program module also 
contains subroutines to perform the tasks selected from the program menu. The last 
task performed by the main program module is to transfer control of the program to 
the clock module. 

The programming for some of the housekeeping and administrative tasks found 
in the main program module are taken in whole or in part from the program written by 
Ε.Ε. Nelsen [Ref. 2]. In particular the subroutine to check for BASICA and the 
graphics card, the subroutine to load the assembly language random number generator 
program, and the subroutine to change the seed of the random number generator were 
taken essentially without change. 

Once a simulation run is initialized and running, the clock module controls the 
flow of the program for the remainder of the simulation run. The first function of the 
clock module is to generate an initial exogenous customer arrival time and determine if 
that customer will exit after completing service at service centers | and 2. This is a one 
time initialization. 

The clock module determines the next-event time as being the smallest of the 
next arrival time, next service completion time at service center 2, or next service 
completion time at service center 2. The next-event time thus determined is always 
rounded up to an integer value so to be compatible with the fixed increment simulation 
clock. 

Every time the simulation clock advances one time increment, the clock module 
updates the state-versus-time graphs for those service centers selected to be displayed 
by the user. Recall that the state of a particular service center is the total number 
waiting and being served by that service center at a particular point in time. When the 
elapsed time reaches 32 time units, the plotted lines on these graphs will be at the right 
most edge of the graphs. When this condition occurs, the clock module shifts the 
displayed state graphs one time unit to the left, and then plots the new time increment. 
This shifting of the graphs by one time unit to the left continues for all successive time 
units. In this manner the service center state-versus-time graphs always contain the 
most recent 32 time units of information, while more dated information is permanently 
lost. When the simulation clock eventually equals the next-event time, the clock 
module shifts control to either the arrival module, service center one module, or service 


center 2 module as appropriate. 
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Every time the clock module increments the simulation clock, it also checks the 
keyboard for certain commands issued by the user. By striking the appropriate 
keyboard keys, the user can suppress the audible tones associated with movement in 
the network, pause the simulation, continue the simulation, change the display 
configuration, speed up the simulation, and end the simulation. These user commands 
are described in detail in Chapter IV, the user's manual. 

There are only five possible events that can happen in the simulation: 

l. An exogenous customer can enter the network. 
2. A customer can complete service at service center | and exit the network. 


A customer can complete service at service center 1 and proceed to 
service center 2. 


4. A customer can complete service at service center 2 and exit the network. 


A Πο can complete service at service center 2 and loop back to service 
center |. 


Event 1 is handled by the arrival module. Events 2 and 3 are handled by the service 
center | module. Events 4 and 5 are handled by the service center 2 module. 

The arrival module is called when an exogenous customer enters the network. It 
flashes the arrow leading into service center 1. If no one is presently in service, a blue 
ball representing the customer is placed in the square box corresponding to service 
center l. If there is a customer already in service but the waiting queue has less than 
eight people in it, a blue ball 1s placed at the end of the queue and the queue counter, 
Q,» is incremented by one. Otherwise only Q, is incremented. 

If the exogenous customer immediately receives service, rather than entering the 
queue, the arrival module calls the random number generator and computes the length 
of service. In all cases the arrival module calls the random number generator and 
computes the inter-arrival time of the next exogenous customer. The inter-arrival time, 
and service time if applicable, are added to the current time (continuous not rounded 
time), to determine when the service will be completed and when the next exogenous 
arrival will take place. When the program returns to the clock module, these newly 
created event times will be compared to the service completion time at service center 2 
to determine what event should occur next in the simulation. 

The service center 1 module is called when the next scheduled event is a service 
completion at service center 1. The departing customer can either depart the network 
Or continue on to service center 2. In either case, the module erases the departing ball, 


shifts any customers waiting for service center 1 to the right, and decrements Q, by one 
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as appropriate. If there is now a customer receiving service at service center 1, the 
random number generator is called and a new service completion time 1s generated. 

If the customer departs the system, the appropriate arrow is flashed. If the 
customer proceeds to service center 2 the arrow leading from service center l to service 
center 2 is flashed. If no one is presently in service at service center 2, the transferred 
ball will be placed in the square box representing service center 2. If there is a 
customer already in service, but the waiting queue in front of service center 2 has less 
than five people in it, the ball 1s placed at the end of the queue, and the queue counter 
Q,, 1s incremented by one. Otherwise, only Q, will be incremented by one. 

If the transferred customer immediately receives service at service center 2, rather 
than exiting or entering the queue in front of service center 2, the service center 2 
module calls the random number generator and determines a service completion time 
for that customer. The service completion time is added to the current time to 
determine when the service will be completed. When the program returns to the clock 
module, this newly created event-time will be compared to the next exogenous arrival 
time and service completion time at service center 2, to determine what event should 
occur next in the simulation. Before returning to the clock module, the service center | 
module calls the random number generator and determines whether or not the next 
customer leaving service center 1 exits the network or proceeds on to service center 2. 

[he service center 2 module is called when the next scheduled event, as 
determined by the clock module, is a service completion at service center 2. The 
departing customer can either exit the network or loop back to service center 1. In 
either case the service center 2 module erases the departing ball, shifts any customers 
waiting for service center 2 to the right, and decrements Q, as appropriate. If there is 
now a customer receiving service at service center 2, the random number generator is 
called and a new service completion time is generated. 

If the customer departs the system, the service center 2 module flashes the 
appropriate arrow. If the customer loops back to service center 1, the bottom arrow 
from service center two to service center one is flashed. If no customer is presently in 
service at service center 1, the transferred ball is placed in the square box representing 
service center |. If there is a customer already in service, but the waiting queue has 
less than 8 people in it, the transferred ball is placed at the end of the queue, and the 


queue counter, Q,, is incremented by one. Otherwise, only Q, is incremented. 
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If the transfered customer immediately receives service at service center |, rather 
than exiting or entering the queue in front of service center 1, the service center 2 
module calls the random number generator and determines a service completion time 
for that customer. The service completion time is added to the current time to 
determine when the service will be completed. When the program returns to the clock 
module, this newly created event time will be compared to the next exogenous arrival 
time and service completion time at service center 1, to determine what event should 
occur next. Before returning to the clock module, the service center 2 module calls the 
random number generator and determines whether or not the next customer leaving 
service center 2 exits the network or loops back to service center Ἱ. 

The service center 2 module performs one other important function. It keeps 
track of how many times a customer loops back to service center 1, allowing this 
number to be physically placed inside the balls representing customers. When a 
customer first arrives and has had no opportunity to loop from service center 2 back to 
service center I, he is depicted as a solid blue ball with no number printed on it. Ifa 
customer has looped from service center 2 to service center | exactly once, he will be 


” 1” 


depicted as a purple ball with a number punted mthe center- ofita Iia customer 
has looped exactly twice, he will be depicted as a purple ball with a number ”2” printed 
in the center of it, etc. 

All of the event modules create audible sounds that correspond to customer 
motion in the network. There are only two distinct tones. A short high pitched tone 
to indicate a departure from the network, and a longer low pitched tone to indicate 
arrivals from both outside the network and arrivals from another service center. 

When a ball arrives at a service center and finds the visible queue space at 
Edipacity ('ò for service center 1“and 5 for service center 2), the modules store these 
customers in a queue that is not visible to the user. These invisible queues, one for each 
service center, preserve the order and identity of the overflow customers. They are 
overflow only in the sense that there is not enough room on the computer screen to 
graphically display them all. If the queue counters (Q,, Q,) do not agree with what is 
physically shown on the screen, it is because there are customers in the invisible queue 
anyqueues, The queue counters are correct. 

After returning to the clock module, the simulation updates the estimates of the 
tabulated performance values, using information generated in the event modules, and 


prints them to the screen. As previously noted, the clock module also updates the 
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service center state-versus-time graphs. Each of the event modules records the 
information necessary for the clock module to know where vertically on the graphs it 
should draw the horizontal line segments representing one time unit. These graphs, 
because of space constraints, can not plot state variables greater than nine. [f the state 
variable is equal to 9, the clock module plots a purple tick mark at a vertical position 
corresponding to 9. For state variables greater than 9, the clock module plots a blue 
tick mark also at a vertical position corresponding to 9. Therefore, observing a blue 


tick mark only indicates that the service center state is 10 or greater. 


G. PROGRAMMING CONSIDERATIONS 

At all times during the simulation, JACKQUE maintains a next-event time for 
each of the three possible events (exogenous arrivals, and service completion at service 
centers 1 or 2). The clock module than takes the minimum of these three times to 
determine the next event. Because of the structure of the network, there are times 
when an event should never be allowed to occur. 

For example, if no customer is currently being served at a service center, the 
program must be prevented from selecting the service completion time of that service 
center as the minimum of the next-event times. If the program did select the service 
completion time as the next-event, the simulation would attempt to move a customer 
from the service center when in fact there is no customer there. 

To prevent the selection of an impossible event, the program sets the next-event 
time of all impossible events equal to 107", which for all practical purposes is equal to 
infinity. The clock module takes the minimum of the next-event times. It will never 
select 1051, because the next exogenous arrival time will always be less than this. In 
this particular model, one or both of the next service times may periodically be set to 
infinity. 

When the clock module selects the minimum of the three next-event times, there 
Is a remote possibility that there will be a two-way, or even three-way tie. In the 
unlikely event that a tie does occur, the program decides randomly which of the events 
should occur first. 

Because of the discrete nature of the graphics display, events can only be 
depicted as having occurred at the end of a discrete time tick. For this reason, within 
one time tick, At, many events may actually occur. Some of these events may occur at 


the very begining of At, some in in the middle, and perhaps some near the end of At. 
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Although their chronological order of occurrence is preserved, they are all plotted 
sequentially as if they occurred at the end of At. For this reason what is observed on 
the screen 1s not an absolutely accurate representation of a Jackson network with 
exponentially distributed inter-arrival times and service times. What 1s graphically 
displayed is a Jackson network with geometrically distributed inter-arrival and service 
times, the geometric distribution being the discrete analog of the continuous 
exponential distribution. Despite this fact, the simulation still produces a realistic 
visual impression of a true Jackson network. 

For a discrete event stochastic process, the state of the process can only change 
by discrete jumps of plus or minus one. However, because JACKQUE considers all 
events occurring in At to have occurred at the end of At, the state-versus-time graphs 
might exhibit jumps greater than plus or minus one. For example if service center | 
experiences 3 arrivals and ] service completion, all during At, the state-versus time 
graph for service center 1 wil jump 2 units when the horizontal line segment 
corresponding to At 1s plotted (see Figure 3.3). Similarly, the state-versus-time graph 
might plot the next tick mark at the same vertical height as the previously plotted tick 
mark, indicating no change when in fact there had been a change during At. For 
example, if service center | experiences one arrival and one service completion during 
At, there is no net change and the time tick will be plotted at the same vertical height 
as before (see Figure 3.4). 

This undesireable effect can be minimized by reducing the probability of having 
multiple events occur during At. This can be accomplished by either reducing the size 
of At or through the judicious selection of the network parameters, as discussed in 
Becton 1V-I. 


H. GENERATION AND USE OF RANDOM NUMBERS. 

As demonstrated by Nelsen [Ref. 2], the random number generator residing in 
BASICA produces pseudo-uniform (0,1) numbers with significant lag-1 correlation. 
This makes it a poor random number generator and unsuitable for use in simulation. 

JACKQUE uses a random number generator modified by Nelsen, and originally 
programmed by B.O. Shubert of the Naval Postgraduate School for use on the 
IBM-PC. This random number generator, named RNGEN, is written in 8088 


assembly language. 
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Figure 3.3 An Example of Jumping 2 States in At. 
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Figure 3.4 An Example of Jumping 0 States in At. 
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RNGEN is a Prime Multiplicative Linear Congruential Generator (PMMLCG) 


based on the following recursive formula: 


Xi = 7° x Mod(2**-1), (eqn 3.1) 
where X; is the old seed and X; 4, is the new seed. A uniform (0,1) random number is 
obtained by dividing the new seed by the modulus. This particular PMMLCG is 
considered to exhibit good statistical behavior and is currently used in the 
LLRANDOM package, IBM’s GGL, and IMSL’s GGUBS [Ref. 6: p.226]. 
JACKQUE uses Uniform (0,1) random variables for three different purposes: 
l. To generate exponentially distributed random variables. 


2. To determine whether a customer exits or remains in the network after 
completing service at a service center. 


3. To break ties in the event two or more event times are minimum. 

The inter-arrival times and service times of customers in a Jackson network are 
assumed to be exponentially distributed. To obtain these exponential random 
variables, Uniform (0,1) random variables produced by RNGEN are transformed 
through the inverse probability method. Let U be a uniformly distributed (0,1) random 
variable, and let X be an exponentially distributed random number with parameter À. 


The distribution function for X is: 
Fx) 2 1-e^* x 2 0. (eqn 3.2) 


An exponential random number is generated by setting u = F(x) and solving for x: 


u= I-e Àx (eqn 3.3) 
l-u = eX (eqn 3.4) 
In(1-u) = -Ax (eqn 3.5) 
x — -l/Aln(1-u) (eqn 3.6) 
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x = -1/Aln(u). (eqn 3.7) 


It is possible to replace (1-u) with u because both (1-u) and u have the same Uniform 
(0,1) distribution. 

When a customer completes service at a service center, he can either exit the 
system or go to the next service center. To determine which occurs JACKQUE 
compares a Uniform (0,1) random number to the fixed probability that the customer 
exits after completing service at service center 1 (P;). If U(0,1) S Pi, the customer is 
considered to have exited. Otherwise the customer 15 considered to have proceeded to 
the next service center: 

A similar procedure is used to break a tie in the event two or three event-times 
are simultaneously minimum. In the case of a two-way tie, event A is chosen if 
U(0,1) S 1/2, and event B is chosen otherwise. In the case of a three-way tie, event A 
is chosen if 0 = U(0,1) < 1/3, event B 1s chosen if 1/3 < UOD € 2/3, and even pena 
chosen otherwise. Because all events do not have equal probabilities of occurring, 
using equal probabilities is not totally accurate, but considered sufficient for the 


purposes of this simulation. 


34 


IV. THE JACKSON NETWORK USER'S GUIDE 


Author 
Gary F. Greene 


James D. Esary 
Thesis Advisor 


SS 


- 


G. 
H, 


TABLE OF CONTENTS 


INTRODUCTION ων a enn touni anti e a e A e e 40 
MODEL DESCRIPTION sos T TM 40 
HARDWARE AND SOFTWARE REQUIREMENTS το... 41 
l..Hardware ον νεος εν ο ο. On e On 4] 
2. Software aa iai komesan n ee A ee ie A a A L ΥΠ 4] 
3. Program Files... uy Te AAA e EDI e en 41 
GETTING STARTED ....... .... NID P M 42 
l. Making an Application Diskette ^ c UI 42 
2. Starting the Simulation «TOT TT n 42 
3. litle Screen. 5... so ον τι 43 
PROGRAM MENU νο το ο. ο τι 45 
l. -Displaying the User Instructions ma nT nn 45 
2. Setting the Random Number Generator seed TT 45 
3. Default.Model] Parameters ο το... T DD 46 
4. Changing the Model Parameters e ου EN 46 
5. Ending the Program. .... 2 n TOT DE 46 
THE GRAPHICS DISPLAY -a PTT 47 
l. The Pictorial Network Display 2 000001 Ne 47 
2. The State-Versus-Time Graph e. το S ο TT 50 
3, Network Performance Measures... 2.2. 2 cee ο U 7 UD 50 
KEYBOARD COMMANDS Sacose rae re ο ο SS SDD πα 32 
PROGRAM LIMITATIONS ;:.2-- 0999 99989 2227 0529) 20072 e 53 
GUIDANCE ON PARAMETER SELECTION Gee ee 54 
PROGRAMING MODIFICATIONS e P pn sn UU D 56 


36 


4.1 
4.2 
4.3 
4.4 
4.5 
4.6 
4.7 


LIS EOF FIGURES 


MA MINE ENN E RE het 44 
MN EI νι ο. Onn ον νεο ee eee eee eee eee ee eee ad 
ου ο ο he eels ae ek n m ye ee ee e eth 45 
ο mol Qí Graphics Display with Unstable Parameters ................... 47 
K DuOnssStiteDiagrad - el ah I RESI] n 48 
Ex pu) ispay With '] woe State Diagrams ................... III 45 
ο πο. o Graphics Screen for Unstable Network ....................... 52 


37 


I. 


LISLOL LABRIEES 


SUGGESTED ΟΠΕ 


9 è è è é+ 9 A è «w *9 « « $# 0€ 6 6 6 a a a © @ @ @ o Q a a Q @ $9 « © @ qq o © @ 


38 


A. INTRODUCTION 

JACKQUE is a discrete event stochastic simulation of a Jackson network. A 
Jackson network is a network of queues where a customer may receive service from 
numerous service centers, possibly visiting the same service center more than once, 
before exiting the network. The program produces a variety of visual displays on the 
computer screen intended to provide the user with a better understanding of the 
dynamic nature of the model. 

This manual is intended to assist the user in the actual use and operation of the 
program. In describing the user interaction with the program, the following notation is 
used: 


e When a command is to be entered (typed) exactly as shown, it will be enclosed in 
single quotes (' ' 


° When. a specific key on the keyboard i Sto Depressed, the letter or letters on the 
key will be encased by brackets (< >). 


For example, 'Load "JACKQUE.BAS",r «Enter», means type the command exactly 
as it appears between the single quotes and then press the ENTER key. PC-DOS and 
JACKQUE make no distinction between upper and lower case letters, so either may be 


used. 


B. MODEL DESCRIPTION 

The Jackson network modeled by this simulation consists of two service centers, 
each containing one server. An arriving customer from outside the network must first 
momuo service center 1 where he enters into a First in First Out (FIFO) queue. Upon 
completion of service at service center 1, the customer either exits the network with 
probability P,, or proceeds on to service center 2 with probability 1-P,, where he enters 
another FIFO queue. Upon completion of service at service center 2, the customer 
either exits the network with probability P,, or loops back to service center 2 with 
probability 1-P,, again entering a FIFO queue. In this model, a customer can loop 
from service center 2 to service center | a large but finite number of times. 

There are several assumptions required for this model. Customers from outside 
the network are assumed to arrive at service center | according to a Poisson process 
with rate A. The service times at service center l and service center 2 are assumed 
exponentially distributed with rates B, and p, respectively. When a customer departs 
the network with probability P, or goes to the next service center with probability 1-P, 


he does so instantaneously. These probabilities are assumed independent of the history 
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of the network. lt is assumed that there is unlimited waiting space at both service 
centers, and that if one server is idle he does not help a busy server at the other service 
CENCEN, 

For this model, the state of the system, N(t), is the sum of iN (t) and N,(t). N;(t) 
(17 1,2) is defined as the number of customers waiting and being served at service 


center i at time t. Therefore N(t) = N,(t) + N,(t) = n, where n — 0,1,2,.... 


C. HARDWARE AND SOFTWARE REQUIREMENTS 
l. Hardware 
An IBM Personal Computer (IBM-PC) with at least 128 Kilo-bytes of random 
access memory (RAM), a color/graphics adapter card, a color monitor, and at least 
one 5 1/4 inch floppy disk drive are required to run the simulation. The color monitor 
and color/graphics adapter are essential since the program, JACKQUE , uses color 
graphics. The 128 K of memory is required because JACKQUE, which 1s written in 
BASICA, uses an assembly language sub-routine which is loaded outside of Basic's 
64 Kilo-byte workspace. 
2. Software 
The primary program JACKQUE.BAS, and the assembly language 
sub-routine RNGEN.SRT are provided to the user on the distribution diskette. In 
addition to these programs, the user must provide the Disk Operating System PC-DOS 
(2.1) and Advanced Basic (BASICA) in order to run the simulation. Because 
JACKQUE directly interacts with the Disk Operating System, there can be no 
guarantee that the simulation will run on JBM compatible machines. 
3. Program Files 
The following program files are included on the distribution diskette: 
e JACKQUE.BAS - the main program written in BASICA. 


e RNGEN.SRT - a random number generator written in 8088 assembly language, 
used as a sub-routine to the main program. 


e JACKQUE.BAT - a batch file that loads and starts the simulation. 


e AUTOEXEC.BAT - an autoexec batch file that loads and starts the simulation 
immediately after the system 1s started or restarted. 


The first two files are absolutely necessary to run the simulation. The two batch files 
are two different methods of automatically loading and starting the simulation. The 
batch files are not essential, but enhance user-friendliness. A source code listing of 
RNGEN.SRT 1s contained in Reference 2. 


40 


P GETTING STARTED 
l. Making an Application Diskette 

The application diskette is created using the distribution diskette. After being 
used for tnis purpose, the distribution diskette should be stored in a safe place, and only 
used to create additional application diskettes as needed. 

Although there are many different ways that a user can create an application 
diskette, this manual will describe just one. A self-starting (bootable) diskette 
containing certain operating files, the BASICA command file and the simulation 
program files is considered the most efficient format. This eliminates the requirement 
for a separate diskette containing PC-DOS, and the diskette can be run using only one 
disk drive. 

To make a bootable diskette, follow the instructions in the PC-DOS manual 
and format a blank diskette using the FORMAT command with the /s option. The /s 
option copies only the DOS files that must be on the newly formatted disk. These are 
the COMMAND.COM file, plus several hidden files that are necessary for DOS to 
operate. The Advanced Basic file, called BASICA.COM, also located on the DOS 
diskette can be transferred onto the application diskette using the copy command. The 
copy command is also used to transfer the simulation’s four files from the distribution 
diskette onto the bootable application diskette. The user now has a totally self 
contained, self starting diskette capable of running the simulation. 

à. Starting the Simulation 

lf the PC is turned off, place the bootable diskette into drive A and turn on 
the computer.” After the computer does it's internal checks taking nearly 60 seconds, 
the computer wil immediately execute the DOS commands contained in the 
AUTOEXEC batch file. These DOS commands automatically load BASICA and 
JACKQUE.BAS, and run the simulation. 

lf the PC is already turned on and in the command level of DOS, place the 
bootable diskette in drive A and type JACKQUE'« ENTER». The computer will 
immediately execute the DOS commands contained in the JACKQUE.BAT file and 
start the simulation. Alternately, with the bootable diskette in drive A, the computer 
can be restarted by simultaneously striking « Ctrl» «alt» «del». In this case the 
simulation is started using the AUTOEXEC.BAT file. 


3This section assumes that drive A is the default drive. 
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[f for some reason the computer is in the command level of DOS and it is not 
desired to use the batch files, the simulation can be started by using the following 
sequence of commands: 

e 'BASICA' « RETURN» 


e 'Load "dJJACKQUE.BAS",r «ENTER? where d is the label of the drive where 
JACKQUE.BAS can be found. 


[f the simulation has been terminated from the program menu, and the 
computer is still in the BASICA environment, another simulation run can be initiated 
by simply typing “RUN’ <ENTER>. If you have returned to the command level of 
DOS, use the methods previously described. 

During the execution of the program, JACKQUE attempts to load the 
random number generator sub-routine from the diskette in the A disk drive. If the 
program fails to find the necessary file in drive A, it asks the user for the correct drive 
label where it can be found. Simply respond by entering the proper drive label and 
<ENTER>. If it can not find the sub-routine on the specified drive, or if the 
specified drive is not valid, the program will again ask the user for the correct drive 
label. This situation will never arise if the bootable application diskette is used in 
drive A. 

3. Title Screen 

Upon starting the simulation, the first thing the user will observe on the 

screen will be the title screen as shown in Figure 4.1. From this point on, the 


on-screen prompts will guide the user through the running of the simulation. 


E. PROGRAM MENU 

Upon leaving the title screen, there will be a few second delay while JACKQUE 
initializes variables and loads the random number generator sub-routine. The next 
screen to be displayed will be the program menu as shown in Figure 4.2. The options 
on the menu allow the user to read the instructions, select a new starting seed for the 
random number generator, run the simulation using the default parameters, run the 
simulation using parameters selected by the user, and end the program. 

When a simulation is ended, JACKQUE always returns to the program menu. 
This allows the user to conduct numerous simulations runs, perhaps with different 
seeds and/or parameters values, during one running of JACKQUE. 

Upon departing the menu, and before running the first simulation, JACKQUE 


quickly plots and then erases the screen shown in Figure 4.3. Before being erased, 
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Press any key to continue... 


Figure 4.1 The Title Screen. 





<I>nstructions 


«N?ew random number Generator Seed 
«D»efault Model Parameters 
<Cohange Model Parameters 

«E»nd the Program 


Enter your selection: 


Figure 4.2. The Program Menu. 


these symbols are stored in memory and are used by JACKQUE in generating graphics 


for all future simulation runs. 


43 


5 
$ 
δ 
e 
3 
dy 
(3 
— 
— 


' Min, 2 
d. ee 














N. ` ` 
Cu m am. 
^ 4 4 4 
` 
a e 
——— - — —P— £ e 
A a 
= ` ` 
A £ “ 
ν ` ` 


Figure 4.3 Graphic Svmbols. 


l. Displaying the User Instructions 

To display the user instructions on the screen, press <I>. The instructions 

displayed are a condensed version of what is contained in this guide. 
2. Setting the Random Number Generator Seed 

A new seed can be selected by pressing <N>. The program will prompt you 
to enter the value of the seed. The seed is entered by selecting a number in the range 
from 1 to 2,147,483.646 and pressing <ENTER>. If anon-integer number is entered, 
it will be truncated to an integer. 

If no seed is selected, a default starting seed is used. If the simulation is 
stopped, and then a new simulation started without leaving JACKQUE and without 
setting a new sced, the generator will continue as if the old simulation had never 
stopped. That is, the first random number generated in the new simulation is the same 
random number that would have been generated next in the old simulation. 

By using the same seed with the same model parameters, the simulation can 
be run numerous times, obtaining the exact same results each time. This may be useful 


for replicating a particular model behavior. 
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3. Default Model Parameters 
If this option Is selected by pressing < D> , the simulation immediately starts 


with the model parameters set to the following default values: 


e À = .l0 
ο" .20 
ο = .l5 
` OU 
ai. —.)5 


If a new seed is desired it must be entered before striking < D>. 
4. Changing the Model Parameters 

When <C> is pressed, the user is prompted for values of the arrival rate, 
service rate at service center 1, service rate at service center 2, probability a customer 
exits after departing service center 1, and probability a customer exits after departing 
service center 2. The rates can be any positive, non-zero real number. The 
probabilities must be between zero and one. The quantities are entered by typing in 
the appropriate value and then pressing <ENTER> in response to each of the five 
questions. If <ENTER> is struck without typing in a number, the parameter value 
will be considered to be Zero. 

For a Jackson network to be stable, the effective arrival rate at each service 
center, taking into account arrivals from both outside the network and from other 
service centers, must be less than the corresponding service rates. If this is not the case 
the network will become unstable in that more customers are joining the queue than 
departing it. 

In the case of an unstable queue, there are no steady state probabilities and 
the network performance measures have no meaning. The stability of the network is 
πας 


simultaneously equal to 0 will always result in an unstable network. Section III-I 


controlled by the interactions of the five model parameters. Setting P 


discusses the selection of model parameters. If the selection of parameters results in an 
unstable network, the screen shown in Figure 4.4 will be displayed. 
The simulation starts immediately after entering the last requested parameter 
value. If a new seed is desired, it must be entered before striking <C>. 
5. Ending the Program 
To end the program, press <E>. This will clear the screen and leave the PC 


in the BASIC command mode. It also returns the screen to the standard monochrome 
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ameter selection results 
ο in s te CEXPLOSIVE) queue. 
The graphical present NO of the 
network will still be accurate, 


r the network statistics 
Will be INULAID, and will NOT be 
displayed. See the user’s manual 


for an explanation. 


C»ontinue with Simulation 
oR Sen to the Program MENU 


ENTER your Selection 





Figure 4.4 Example of Unstable Parameter Warning. 


mode with column length equal to 80. If it 1s desired to return to the DOS command 
level type system < ESTERE 


F. THE GRAPHICS DISLPEAY 

JACKQUE uses various graphical displays to allow the user to observe the flow 
of customers through the network. When the simulation is first started, the graphics 
display will appear as shown in Figure 4.5. The graphics display is divided into three 
regions: a pictorial display of the network, the state-versus-time graph for service 
center I, and a table of network performance measures. 

At any time after the start of the simulation, the user can chose to change the 
configuration of the graphics display to that shown in Figure 4.6. The table of 
performance measures has been replaced by the state-versus-time graph for service 
center 2. The user can freely change the configuration back and forth between the two 
modes as described in Section IV-G. The following descriptions are directed at the 
display mode shown in Figure 4.5. 

l. The Pictorial Network Display 

The display in the upper half of the screen is a pictorial representation of the 
Jackson network at a particular instant in time. The left square box represents service 


center l. The right square box represents service center 2. The vertical and horizontal 
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Figure 4.5 Screen Display with One State Diagram. 
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Figure 4.6 Screen Display with Two State Diagrams. 


lines connecting these two boxes represent the possible paths that a customer can 


follow through the network. An arrow directed into a service center is the path of an 
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arriving customer. An arrow directed out of a service center whose head is 
unconnected, is the path of a customer departing the network. Arrows connecting the 
two service centers are the paths of customers moving from one service center to the 
other. The lowest horizontal line directed from right to left is the feedback loop from 
service center 2 1ο ο Ίος Gen mls 

Customers are represented as colored solid balls. A customer arriving from 
outside the network for the first time is represented as a solid blue ball. If a customer 
loops from service center 2 to service center 1l, the color of the ball is changed to 
purple, and a number "1^ is placed inside the ball. Each time this customer loops back 
from service center 2 to service center 1 the number inside the purple ball is 
incremented by one, up to a maximum of nine. Thus a purple ball with a number ^57 
in it, represents a customer who has been fed back from service center 2 to service 
center | five times. When a customer exits the network, the corresponding ball 
disappears from the screen. 

A customer actually receiving service from a service center 1s placed inside the 
appropriate square box. If a customer is waiting for service at a service center, he is 
placed in the waiting line immediately to the left of the appropriate service center 
(square box). The counters under the two waiting lines, Q, and Q,, display the 
number of customers waiting for the respective service centers. These counters do not 
include the person actually receiving service. Because of the limited size of the 
computer screen, it is possible for the values of Q, and Q, to be greater than the 
number of customers actually shown waiting as balls. These invisible waiting customers 
are still accounted for, and will enter the visible waiting line when space permits. 

The five model parameters, À, jt,, jt,, P,, and P, are displayed in appropriate 
places on the pictorial network display. The current time of the system is printed at 
the bottom of this display. 

Each time a customer (ball) moves in the network, the motion is represented 
by moving the affected balls, flashing the appropriate path arrow, and sounding a 
distinctive tone. There are only two distinctive tones. A short high pitched tone to 
indicate a departure from the network, and a long low pitched tone to indicate arrivals. 
Arrivals consists of both arrivals from outside the network and arrivals from other 
service centers. 

In Figure 4.5, both service centers are busy with four customers waiting for 


service center 2 (Q,=4), and twelve customers waiting for service centers o ey 
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although only eight are physically shown. Of the visible customers, one has 
experienced feedback twice, four have experienced feedback once, and nine are passing 
through the network for the first time. 

2. The State-Versus- Time Graph 

The lower left corner of the graphics display contains the state-versus-time 
graph for service center 1. The values of the state variable N(t), the number of 
customers waiting and being served at service center | at time t, are plotted on the 
ordinate. Time is plotted on the abscissa, and is represented as 32 discrete divisions. 
Each of these time segments corresponds to one time tick of the simulation clock. 

After each increment of the simulation clock, the program determines the 
current state of service center l. The program then uses this information to plot a 
horizontal line one time tick in length at the appropriate height on the 
state-versus-time graph. In Figure 4.5, at time 3818, the state of service center 1 went 
from 8 to 7. 

When the elapsed time reaches 32 time units, the plotted line on this graph 
will be at the right most edge of the graph. When this condition occurs, the program 
shifts the entire graph one time unit to the left, and plots the new time increment at the 
appropriate height. In this manner the graph always contains the most recent 32 time 
units of information on the state of service center 1. An identical description pertains 
to the state-versus-time graph for service center 2, which can be displayed at the users 
option. 

The state-versus-time graphs, because of space constraints, can not plot state 
variables greater than 9. If the state variable is equal to 9, a purple tick mark at a 
vertical position corresponding to 9 is plotted. For state variables with values greater 
than nine, a blue tick mark 1s plotted, also at a vertical position corresponding to 9. 
Therefore observing a blue tick mark only indicates that the service center state is 10 or 
greater. 

3. Network Performance Measures 

The lower right hand corner of the graphics display contains a table of 
estimated and theoretical equilibrium quantities. These quantities are just a few of the 
many that can be used to evaluate the performance characteristics of the network. The 


variable names used in this table are defined as: 


R; = the effective arrival rate of customers at service center 1 (17 1,2). 
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S. — the expected number of times a customer will pass through service 
l center i before exiting the network (i= 1,2). 


L; = the expected number of customers both waiting and being 
Served at service centes rsi 


L = the expected number of customers in the network, both waiting 
and being served. 


W. = the expected amount of time a customer spends waiting and being 
served at service center 1 (1= 1,2). 


W = the total expected amount of time a customer spends in the network, both 
waiting and being served. 


Z = the proportion of time service center | is busy and service 
center 2 is idle: 

The quantities that appear under the column header Limit, are the limiting 
values of these quantities. They are computed using results from Little's Law and from 
the analysis of simple M/M/1 queues. They are the theoretical values that would be 
obtained if the network reached steady state. 

The quantities that appear under the column header Estimate are estimates of 
these theoretical limiting values. They are computed with data collected from the 
simulation. These estimates are updated after every event occurrence, where events are 
defined to be the arrival of a customer from outside the network and service 
completions at either service center. 

At certain points during the simulation, some of the estimates might be 
negative and/or extremely small or large in magnitude. If an estimate exceeds the 
format capabilities of JACKQUE by being too large or small, the estimate will be 
displayed on the screen as "9999,99", Some of the estimates will temporarily become 
negative if the estimated effective arrival rate at a service center 15 greater than the 
estimated service rate of the service center. To avoid confusion, the absolute value of a 
negative estimate will be displayed on the screen. 


Because of the manner in which Lol, Wo ane Z are estimated, the 


2^ 
estimates might show considerable variance around the theoretical value. This 
variability will decrease as the simulation progresses in time. 

lf in the selection of model parameters, an effective arrival rate 1s greater than 
the service rate at the applicable service center, the user 1s giyen OTa 


terminate the current simulation or to continue as described in Section IV-E-4. If the 
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user opts to continue with the present model parameters, the resulting theoretical and 


estimated values will no longer be valid. In this case the graphics display will appear 


as seen m Figure 4.7. 
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Figure 4.7. Example of Graphics Screen for Unstable Network. 


KEYBOARD COMMANDS 


While the simulation is running, the user has the ability to interact with the 


model through the use of the following keyboard commands: 


« P» - pressing P for “Pause” will temporarily freeze the simulation, allowing the 
user time to study the graphical display. 


<C> - pressing C for “Continue” counteracts pressing <P> and resumes the 
simulation where it left off. 


<B> - pressing B for “Beep” will turn_off the audible tones that accompany 
customer motion through the network. Pressing <B> a second time will turn 
the beeps back on. 


<D> - pressing D for “Display” replaces the table of estimated quantities with 
the state-versus-time diagram for service center 2. Pressing <D> a second time 
returns the graphics screen to the original configuration. 


<F> - pressing F for “Fast”, accelerates the simulation, causing the estimated 
values to sone to the theoretical values more πο In the fast mode, the 
pictorial network display and state-versus-time graph for service center | are 
erased and not updated. Additionally the estimates are only updated after every 
10 event occurrences. Pressing < F> a second time slows the simulation and 
resumes ο ms the graphics. Upon return to the slow mode, it may take some 
time τς the pictorial network display to catch up with the actual status of the 
network. 
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e «S» - pressing S for “Stop”, returns the simulation to the program menu. 
Because JACKQUE only looks for keyboard commands while incrementing the 
simulation clock, there may be a slight delay between the user input and the program 


response. 


H. PROGRAM LIMITATIONS 

JACKQUE always starts the simulation with the network completely empty of 
customers. There is no way of starting the simulation in any other state. It 1s also not 
possible to change the model parameters once the simulation is started. If it is desired 
to do so, the present simulation must be terminated, the parameters changed, and a 
new simulation started. 

JACKQUE uses string arrays to keep track of the customers and their feedback 
count as they move through the network. The maximum length of a string in BASIC 
is 256. Because of this restriction, JACKQUE can not have more than 256 customers 
Waiting for either service center | or service center 2. If this condition is met, the 
simulation will automatically be terminated and the following message displayed on the 
Screen: 

Program execution HALTED. Queue 
length exceeds program capabilities. 
Press any key to return to Main 
Program Menu. 

JACKQUE 1s incapable of explicitly tracking any individual customer for more 
than 9 feedbacks from service center 2 to service center l. Any customer who has 
experienced feedback nine or more times will be represented as a purple ball with à 
number ^9" in the middle. Upon seeing this symbol, it is not known if that customer 
has experienced feedback 9 or more times. This limitation causes the estimates of S, 
and S, to be smaller than the theoretical values. This bias can be minimized through 
the judicious selection of the model parameters as discussed in Section IV-I. 

Although model events are occurring continuously in time, JACKQUE updates 
the graphics display only after the passage of a discrete amount of time, At. It is 
therefore possible for one or more events to occur during At. When JACKQUE plots 
the events on the pictorial display, all events occurring during At are plotted as if they 
occurred at the end of At. They are plotted one immediately after the other, in the 
proper chronological order in which they actually occurred. 

The time it takes the computer to physically plot the movement of the balls is 


not considered by JACKQUE as part of the simulation time. Because the number of 
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events plotted after each At is variable, and because certain events take longer to plot 
than others, the illusion of real time passage to the user is some what distorted. By 
distorted, it is meant that the increments of the clock as shown on the graphics display 
are not always of equal time duration. 

The illusion of real time passage can be maintained, and the frequency of 
multiple state jumps can be minimized through the judicious selection of the model 


parameters as discussed in the next section. 


ll GUIDANCE ON PARAMETER SELECTION 

The behavior of the network is controlled by the five parameters À, μ., p 
and P,. Many of the shortcomings of the graphical display can be overcome by the 
careful selection of these parameters. 

The illusion of real time passage can be preserved by limiting the frequency with 
which multiple events occur during At. The inter-arrival rate and service rates should 
be selected small enough so that the probability of numerous arrivals or service 
completions occurring in At is small. Picking small rates will also improve the users 
ability to follow the progress of a customer in the network. If the customers are 
moving too quickly, the graphics display can become confusing to the user. Finally, 
the selection of small rates will minimize the number of times that the state-versus-time 
graphs will exhibit multiple jumps. 

Although some feedback makes the model more interesting, it 1s desirable to keep 
the average amount of looping for a customer below nine. This minimizes the amount 
of bias introduced into the estimates for S, and S,. It also reduces the number of 
times a customer is represented as a ball with a number “9” in the center. This is 
desirable because the user does not know the true meaning of a number nine. The 
amount of looping is completely controlled by the values selected for P, and P,. 

It is also desirable to keep the average queue size in front of the service centers 
close to the number that can physically be shown on the computer screen. Although 
this is not essential, it makes the graphic simulation easier to understand. This is 
accomplished by keeping the service rates at the service centers some what larger than 
the effective arrival rates. 

Selecting the parameters to accomplish all of these things can be difficult because 
they are so highly interrelated. Table 1 contains some parameter values that have 
yielded interesting results. One possible method for selecting parameters is presented 


below: 
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51 
It is 


Select À, the arrival rate of customers arriving at service center | from outside 
the network. 


Select P, and P,, the probabilities of exiting the network at service completion. 
Compute R, and R,, the effective arrival rates: 

R. λ/{1-{1-Ρ.Χ1-Ρ.}) 

R, ΠΠ.) 


Select jt, and ji, so that: 
H, > F, 


li, R, 


Run the simulation with these parameters and see how it performs. 


required that p, and p, be greater than zero so that when calculating the 


theoretical performance quantities, division by zero 1s not attempted. 





TABLE I 
SUGGESTED PARAMETERS 


Pa Hy μ; Š, 
0.06 0.05 55 
0.14 0.12 2.8 
0.60 0.50 5) 
0.25 0.20 2:3 
0.60 0.55 59 
0.15 0.07 lè 
0.16 0.05 1.0 
0.50 0.35 2:3 
0.60 0.36 ES 


The Jackson network being modeled can be transformed into a simple M/M/1 


queue by selecting P, = 1. Likewise, it can be transformed into a simple M/M/2 queue 


by selecting P, = 0 and πα. 
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J. PROGRAMING MODIFICATIONS 

This section is intended to describe some minor changes that can be made to the 
program JACKQUE.BAS. A BASIC programing manual should be consulted for the 
specifics of how to change and then save the altered program. 

Program line 3306 specifies the default drive on which the program looks for the 
random number generator sub-routine RNGEN.SRT. To change the default drive, 
simply replace DRIVES = ^A:" with the appropriate drive label. For example, if it is 
replaced with DRIVES="B:”, the program will automatically look on drive B for the 
sub-routine. 

Program line 3455 contains the default parameters. These are the parameters 
used in the simulation if the user does not specify something else. To change one or 
more of the default parameters, simply replace the existing value with the desired value. 

The lengths of the audible tones used in the program are set in program line 
3456. The variable DURI sets the duration of the low pitched tone corresponding to a 
customer arrival. The variable DUR2 sets the duration of the high pitched tone 
corresponding to a departure from the network. Either tone duration can be changed 
by simply assigning a different positive integer to the appropriate variable. The 
variable assignment must also be changed in program lines 8605 and 8610. 

The actual tones of the sounds are determined by the number appearing in the 
statements “SOUND 3000, DUR2”, and “SOUND 2500,DUR1”. To change the tones, 
simply change the number appearing after SOUND. These statements appear in 
program lines 12035, 18015, 18135, 23015, and 23135. 

The relationship W 2 N in program line 5088 determines how often the 
estimated quantities in the displayed table are updated. They are updated every time 
W exceeds N. W is the count of how many events have occurred since the last update. 
The value of N is set in lines 3470, 9505, and 9515. The variable W is initialized in 
program line statements 3475 and 5088 and incremented in line 5156. As presently 
Written, in the slow mode the estimates are updated after every event and in the fast 


mode they are updated after every 10 events. 
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V. CONCLUSIONS 


This paper presented the development of a graphic simulation of a Jackson 


network. The method of graphic simulation allows the user insight into the dynamic 


nature of the Jackson network, not possible through other instructional methods. 


There are numerous other systems that could be better understood through the 


use of graphic simulation. A few that come to mind are: 


Queues with finite storage. 

Queues with bulk arrivals 

Queues with bulk service. 

Multi-channel queues with one waiting line. 

Multi-channel queues with multiple waiting lines, with or without jockeying. 
Priority queuing. 

Queues with balking customers. 


It is hoped that this effort has been able to assist the reader in better 


understanding Jackson networks. 
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1000 
1005 
1010 
1015 
1020 
1025 
1030 
1035 
1040 
1045 
1046 
1050 
1051 
2000 
2005 
2010 
2015 
2020 
2025 
2030 
2035 
2040 
2045 
2050 
2055 
2060 
2061 
2100 
2111 
2112 
2113 
2114 
2120 
2130 


2140 


2150 


2160 


APPENDIX A 
JACKQUE.BAS PROGRAMI LISTING 


'++++T+TiTliTliTlYT1T17vV7171T1++T+++T++T++T++++++++i+iflitit+Tii+++++T+1T++T++7+7+7++++++++++++++++++++itti++++ 
1+++++++++++Ti+++ 
abb BR BR GÀ B RR  À . À B À 1| 4 4 4 4 B 4 6 À À BG BÀÁERBRG B Y À Á BÀ B BG BGAB GA BO B BG B GÀ BG BG BG BG GÀ BA BGARURRER BEBE BGÀ BA À M G GÀ GB DB AIABA4AA 
'JACKQUE.BAS, Ver 1.0, G. F. Greene; 
CLEAR,,7300: KEY OFF: CLS 

GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOTO 5000 
SCREEN 0: WIDTH 80: PRINT "Program Ended at Your Request..." 
END 


MAIN PROGRAM MODULE Z MPM .++++++++f++i+++i++i+++++++++++i+i++++ 


Aug 1986 


1000 
2000 
2100 
3200 
3300 
3400 
3500 
3600 
3900 


‘check for basica and color graphics card 
‘Title screen 

'One time initialization and Housekeeping 
'load user defined functions 

‘load random number generator 

‘initialization for individual simulation run 
‘print main program menu 

‘compute theoretical statistics 

‘print main display 

‘transfer control to Clock Module 


======= BASICA & COLOR/GRAPHICS ADAPTER CHECK SBRT / MPM 
DEF SEG = O: IF (PEEK(8H410) AND &H30) <> &H30 THEN DEF SEG: GOTO 1030 
LOCATE 3,1: PRINT "Sorry...." 
PRINT "You do not have the color/graphics monitor adapter!" 
PRINT “This simulation uses graphics and requires that adapter." 
DEF SEG: END 
ON ERROR GOTO 1035: PLAY "pló": GOTO 1050 
WIDTH 80: CLS: LOCATE 3,1 
PRINT “This simulation uses advanced BASIC." 
PRINT "Return to DOS by typing SYSTEM‘ and then" 
PRINT "reload this program after using the command BASICA'.":END 
ON ERROR GOTO 0: RETURN 'return to line 60 


TITLE SCREEN / MPM 


SCREEN 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


SCREEN 


0,1: COLOR 14,1,1: WIDTH 40: CLS 


2212: PRINT 
4,19: PRINT 
6,12: PRINT 
9,16: PRINT 
14,5: PRINT 
15,2: PRINT 
17,1: PRINT 
19,4: PRINT 
20,11: 
2224: PRINT 
25,2: PRINT 


= ONE TIME INITIALIZATION AND HOUSEKEEPING SBRT / MPM 


"GRAPHIC SIMULATION" 

"of a" 

"JACKSON QUEUE NETWORK" 

"by G. F. Greene" 

"Submitted in partial fulfillment" 

"Of the requirements for the degree of" 
"Master of Science in Operations Research" 
"from the Naval Postgraduate School" 


PRINT "Monterey, California" 


"Advisors J.D. Esary» G.F. Lindsay" 
"Press any key to continue..."5;: GOSUB 29000:RETURN 


0,1: COLOR 14,1,1: WIDTH 40: CLS 


LOCATE 
LOCATE 


10,5: PRINT "Program Initialization . . . " 

12,13: PRINT "Please Wait" 

LOCATE 13,13,0 ‘hide the curser 

DEFINT A-Z ‘define all variables as integers 

DIM STAT!(2), Q( 2), NSCT!(2), EXIT(2), STRGI$(255), STRG2$(255),SERV2(2200); 
CIRCO(35), CIRC1(35), CIRC2(35), CIRC3(35), CIRC4(35),circ5 (35), CIRC6(35), 
CIRC7(35), CIRC8(35), CIRC9(35),TEMPCIRC(35) 

DIM ARROW1(300), ARROW1C(300), ARROW2(100), ARROW2C( 100), ARROW3(300)>, 
ARROW3C(300), ARRONG( 200), ARROWGC( 200), ARROW5(600), ARROW5SC(600),SHIFT( 2200), 
STATISTICS( 2200), SERVBY1(10), SERVBY2(10) 

FMT1$ = “8888: FMT2$ = “HHH. HHH : FMT3$S = "H.R": 
FMT5$ = "#88" 

TRUE = -1 : FALSE = 


FMTGS = "##.##": 


0: FIRSTSCREEN = TRUE 
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2165 
2166 
3200 
3204 
3205 
3210 


3220 
3221 
3300 
3305 
3306 
3310 
3315 
3320 
3325 
35207 
3329 
3330 
3335 
3336 
3345 
3350 
3355 
3356 
3400 
3420 
3425 
3430 
3440 
3445 
3446 
3447 
3448 
3449 
3450 
3451 
3452 
3453 
3454 
3455 
3456 
3457 
3458 
3459 
3460 
3461 
3462 
3463 
3464 
3465 
3469 
3470 
3475 
3485 
3486 
3500 
3505 
3510 
$515 
3520 
3525 
3535 
3540 
3545 
3547 
3550 
3555 
3560 


RETURN ‘return to line 80 


Iz==z==========s======= USER DEFINED FUNCTIONS SBR 7 NEM. -------------------— 

DEF FNROUNDUPZ(A!) = INT(A!) + 1 

DEF FNEXPONT!(AI!, A2!) = -LOG(A1! )ZAZ! 

DEF FNMIN!(A!,B!,C?) = (-(A!<B! AND A!SCT)XA!) + (-(B!SA! AND BISC! )XB!) + 
(-(C* «A! AND C'«B! )xC!) 

RETURN  'return to line 90 


'zzzz2zzzzzzz2zzczzc2 LOAD RANDOM NUMBER GENERATOR / MPM szzszzzzzzzzzzzzzzzzz- 
RNGEN-0: U!-0: 'initialize random number generator 

DRIVES = "A:" ‘program looks on A drive for RNGEN.SRT first 

ON ERROR GOTO 3329 'RNGEN.SRT cannot be found on default disk drive. 

DEF SEG = &H1A00: BLOAD DRIVES * "rngen.srt",0 ‘load rngen.srt 

ON ERROR GOTO O ‘turn error trapping off 

RETURN ‘return to line 100 


ΙΣΞΞΕΞΕΞΞΞΞΞΞΞΞΞΞΕΤΙΕ LOADING ERROR TRAPPING SBRT ⁄/ MPM ============s=s==s=s= 
IF ERR <> 53 AND ERR <> 76 THEN ON ERROR GOTO O ‘turn error trapping off 
CLS: LOCATE 3,1: PRINT "On which drive can the random number" 

PRINT "generator program be found?" 

GOSUB 29000 ‘capture Keyboard input 

DRIVES = RESPONSES + ":" 

RESUME 3315 


'zzzzz--zz- INITIALIZATION OF SINGLE SIMULATION RUN SBRT / MPM ========== 
FOR I = 1 TO 2 


Q(I) = 0 ‘queue counters 
STAT(I) = 0 ‘status of servers: l=busy O=idle 
NSCT'(I) = 1E+Z6 ‘set next service time at 1 equal to infinity 
NEXT 
FOR I = 1 TO 10 
SERVBY1(I) = 0 ‘number of customers served I times by server 1/2 
SERVBY2(I) = 0 
NEXT 
MODE = 2 ‘graphics mode: 2=with stats table} l=w/o table 
TIME = 0 ‘integer clock time 
CTIME! = 0 ‘continuous real time 
TA] = 0 ‘total arrivals at server 1 
TA2 = 0 ‘total arrivals at server 2 
LAMBDA! = .1 : U1! = .2 : U2! = .13: Pl! =.3 : P2! =.5 ‘default parameters 
DUR1 = 3 : DUR2 = 2 ‘duration of audible beeps 
UNSTABLE = 0 ‘unstable parameters O=no l=yes 
El = O 'total exits after server 1 
E2 = O 'total exits after server 2 
ET = O 'ET = E1 + E2 
NA1 = O 'exogenous arrivals at server 1 
FB1 = O 'feed-back arrivals at server 1 
FAST = 0 ‘display speed: O=slow l=fast 
CUMSERV1! = 0: CUMSERV2! = 0 ‘amount of time server 1 & 2 busy 
STRG1$ = "" ; STRG2$ = "* ‘queue strings for server 1 and 2 
NEXTEVENT = 0: XX = 23 
N = 1 'how frequently stats table updated 
W= 0 ‘event counter 
RETURN ‘return to line 110 
aS === m MAIN PROGRAM MENU SBRT 7/7 MPM zszszszszszszzzzzzzzzzzzzzz 
BORDERS = STRINGS$(40,223 ) 


SCREEN 0,1: COLOR 14,3,0: WIDTH 40: CLS 

LOCATE 3514: PRINT "PROGRAM MENU'': PRINT BORDERS 
LOCATE 6,3: PRINT "<I>nstructions" 

LOCATE 8,3: PRINT "«N»ew random number Generator Seed" 
LOCATE 10,3: PRINT "«D»efault Model Parameters" 
LOCATE 12,3: PRINT "«C»hange Model Parameters" 
LOCATE 14,3: PRINT "«E»nd the Program" 

LOCATE 15,3: PRINT BORDERS 

LOCATE 18,3: PRINT "Enter your selection:"'; 
GOSUB 29000 'capture Keyboard response 

AZ = INSTR("INDCEindce", RESPONSES ) 
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3565 


3566 
3600 
3605 
3610 
3620 
3625 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 
3720 
3721 
3900 
3905 
3910 
3915 
3920 
5955 
3940 
3945 
3950 
3955 
3960 
3965 
3970 
3980 
3981 
4000 
4001 
4002 
4003 
4004 
4006 
4010 
4012 
4014 
4016 
4018 
4020 
4022 
4024 
4025 
4028 
4030 
4032 
4034 
4036 
4038 
4040 
4050 
4051 
4052 
4054 
4056 
4058 
4060 
4062 
4064 
4065 
4066 
4067 
4068 


IF AZ = O THEN 3510 ELSE ON A7 GOTO 4000,4200,120,4500,150,4000,4200 
;120,4500,150 


EEIZIE-CEECEECS THEORETICALSSTATISTICS SBRT / MPM zzzzzzzzzzzz2zz2zzzz-- 
IF (P1'z0 AND P1'!'zP2!') THEN UNSTABLE =1: GOTO 4600 ‘explosive queue 
RI! = LAMBDAt/(1-(1-P1! }*(1-P2!))} 

R2! = (1 - Plt) * RI! 

IF (R1'»-U1'! OR R2!'>=U2!) THEN UNSTABLE= 1: GOTO 4600 ‘explosive queue 
S1! - 1/(P1! * ((1-P1!')*(P2!'))) 

S2! (1-P1!)* S1! 


LI! - R1'/(U1! - R1!) 
L2* - R2'!'/(U2' - R2!) 
L! = L1! + L2! 
W1! = L1! /R1! 
W2' - L2'/R2' 


W! = L!/LAMBDA! 

Z' - (R1'!/U1!)*(1-R2!/U2'!) 
RETURN ‘return to line 130 
' 


'zzzzzzzzzzzzzzz PRINT MAIN DISPLAY SBRT / MPM sszszzzzzzzzzzzzzzzzzzzzz 
SCREEN 1,0: COLOR 0,1 

IF FIRSTSCREEN THEN GOSUB 4700: FIRSTSCREEN = FALSE ‘one time printing 
GOSUB 4800 'main screen drawing sbrt of graphics symbols 
GET (158,93) - (319,199), STATISTICS ‘put stats table on display 
LOCATE 11,21: PRINT USING FMT1$3 TIME 

LOCATE 3,8 : PRINT USING FMT4$3 LAMBDA! 

LOCATE 9,19: PRINT USING FMTG$}3 U1! 

LOCATE 9,33: PRINT USING FMT4$}3 U2! 

LOCATE 3,26: PRINT USING FMT3$3 Pl! 

LOCATE 3,37: PRINT USING FMT3$; P2! 

LOCATE 7,27 : PRINT USING FMT5$; Q(2) 

LOCATE 7,10 : PRINT USING FMT5$; Q(1) 

RETURN 'return to line 140 


!'zzzzzzzzzzzzzzzzzzz USER INSTRUCTION SBRT / MPM ΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞ55Ξ555 
SCREEN 0,1: COLOR 14,1,1: WIDTH 40: CLS 

LOCATE 1,6: PRINT “Instructions for the Program Menu" 

LOCATE 3,2: PRINT "Press <s> to set the random number seed"; 
LOCATE 4,1: PRINT "If «d» or «c» is pressed before setting" 
LOCATE 5,1: PRINT "a seed, the default seed will be used" 
LOCATE 7,2: PRINT "Press «c» to change model parameters." 
LOCATE 8,1: PRINT "The program will ask for the arrival" 
LOCATE 9,1: PRINT "rate (lambda), the service rate for" 
LOCATE 10,;,1:PRINT "server 1 (ul), the service rate for" 
LOCATE 11,1:PRINT "server 2 (u2), the probability of exit" 


LOCATE 12,1:PRINT "after server l (P1), and the probability" 
LOCATE 13,1:PRINT "of exit after server 2. The rates can" 
LOCATE 14,1:PRINT "be any positive non-zero number while " 
LOCATE 15,1:PRINT "P1 & P2 must be between O and 1. The" 
LOCATE 16,1:PRINT "simulation starts after entering P2." 
LOCATE 18,2:PRINT "Press «d» to use the default parameters"; 
LOCATE 19,1:PRINT "The simulation will immediately start'; 
LOCATE 20,1:PRINT "with the values lambda=.10, ul=.20," 
LOCATE 21,1:PRINT “u2=.13, P1=.30, and P2=.50."' 

LOCATE 23,2:PRINT "Press <e> to end the program." 

LOCATE 25,2:PRINT "Press any key to continue..." 3: GOSUB 29000 
CLS: LOCATE 1,13: PRINT "Keyboard Commands" 

LOCATE 3,2: PRINT “During the simulation you may use the" 
LOCATE 4,2: PRINT "following Keys:" 

LOCATE 6,4: PRINT “<p> - Pause the simulation." 

LOCATE 8,6: PRINT "«c» - Continue the simulation." 

LOCATE 10,4:PRINT "«b» - turn on and off the audio Beeps" 
LOCATE 12,4:PRINT "«d» - change the Display" 

LOCATE 13,10:PRINT “configuration back and forth " 

LOCATE 14,10:PRINT “between the state diagrams and" 

LOCATE 15,10:PRINT "the statistics table." 

LOCATE 17,4: PRINT "«f» - turn on and off the" 

LOCATE 18,10: PRINT "fast display mode." 

LOCATE 20,4:PRINT "«s» - Stop the simulation and " 
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4070 
4072 
4080 
4081 
4082 
4086 
4086 
4088 
4090 
4092 
4094 
4096 
4098 
4100 
4102 
4104 
4106 
4120 
4122 
4124 
4126 
4128 
4130 
4132 
4136 
4136 
4138 
4160 
4142 
4144 
4146 
4148 
4150 
4152 
4160 
4162 
4164 
4166 
4168 
4170 
4172 
4176 
4176 
4178 
4180 
4182 
4184 
4186 
4188 
4190 
4192 
4196 
4196 
4198 
4199 
4200 
4205 
4210 
4215 
4220 
4225 
4230 
4235 
4240 
4245 
4250 
4255 
4260 
4265 


LOCATE 21,10:PRINT "return to the program menu." 
LOCATE 25,2:PRINT “Press any Key to continue...'5: GOSUB 29000 


CLS: 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


LOCATE 1,13: 


3,2: PRINT 
4,1: PRINT 
6,1: PRINT 
7,7: PRINT 
9,1: PRINT 
11,1: PRINT 
13,1: PRINT 
15,1: PRINT 
17,1: PRINT 
18 ,5 : PRINT 
20,1: PRINT 
21,5: PRINT 
23»1: PRINT 
25,2:PRINT 


PRINT "Graphics Display" 

"The variables appearing on the display 
"are defined as:" 

"lambda: Arrival rate from outside the" 
"queue network." 

"ul: Service rate for server 1." 

"142: Service rate for server 2." 

"Pl: Probability of exit after server 1." 
"P2: Probability of exit after server 2." 
"Q1: Number of people waiting for" 
"server 1." 

"Q2: Number of people waiting for" 
"server 2." 

"Time: Current time on simulation clock." 
"Press any Key to continue...'$;$: GOSUB 29000 


CLS: LOCATE 1,1: PRINT "N(t): No. of people waiting and/or being 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


CLS: LOCATE 1,1: PRINT "L: L = 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


GOTO 3500 


2,7: PRINT 
4,1: PRINT 
5,5: PRINT 
7,1: PRINT 
8,5: PRINT 
10,1:PRINT 
11,5: PRINT 
13,1: PRINT 
16,5: PRINT 
16 ,1: PRINT 
17,5: PRINT 
18 ,5:PRINT 
20,1: PRINT 
21,5: PRINT 
2255: PRINT 
25,2:PRINT 


3,1: PRINT 
4,5: PRINT 
6,1: PRINT 
7,5: PRINT 
9,1: PRINT 
11,1: PRINT 
12,5: PRINT 
14,1: PRINT 
15,8: PRINT 
16,8: PRINT 
17,8: PRINT 
19,1: PRINT 


"served by given server at time t." 

"Rl: Effective arrival rate at server 1," 
"including feedback." 

"R2: Effective arrival rate at server 2," 
"including feedback." 

"S1: Expected no. of times 
"be serviced by server 1." 
"S2: Expected no. of times a person will" 
"be serviced by server 2." 

"L1: Long run expected no. of people" 
"either waiting or being served by" 

"server 1l." 

"L2: Long run expected no. of people 

"either waiting or being served by" 

"server 2." 

"Press any Key to continue..."3: GOSUB 29000 
Lietz 

"Wl: Expected time a person 
"waiting or being served by 
"W2: Expected time a person 
"waiting or being served by 
"W: Total expected time in network." 

"Z: Long run proportion of time server 1" 
“is busy and server 2 is idle." 

“Limit: Values listed under the column" 
"heading limit' are the limiting" 

“values as time goes to infinity" 

“of the above defined quantities." 
“Estimate: Values appearing under the" 


a person will" 


spends" 
server 1." 
spends" 
server 2." 


20,11:PRINT"heading estimate’ are the" 
21,11:PRINT"estimates of the limiting" 
22,11:PRINT"values. They are calculated" 
23,11:PRINT"with data collected from" 

26 ,11:PRINT"the simulation.",; 


£5, 2:PRINT 


COLOR 195,050: CLS 


RANDOM NUMBER GENERATOR SEED SBRT / MPM 


"Press any Key to return to the menu...";:GOSUB 29000 


'return to the menu 


LOCATE 1,3: PRINT "Set Random Number Generator Seed": PRINT BORDER$ 
LOCATE 4,2: PRINT "Permissible seed values are integers" 

PRINT " in the range: 1 to 2147482646." 

LOCATE 7,2: INPUT "Enter the seed value : " ,A# 

IF A#<1 OR A#>2147483646# THEN LOCATE 7,2: PRINT STRING$(39,32): GOTO 4225 
DEF SEG = &H1A00 

Al! = INT(A#/16777210#): A2! = INT( (A#-A1!¥*167772108% )/65536! ) 

‘poke the seed's upper 2 bytes into RNGEN's seed storage. 

POKE &H164,Al1!: POKE &H163,A2! 

AR=AR-A1'*%167772108#-AZ2!*65536!: AL! = INT(AR/256): AZ! = AR-A1!'*256 
‘poke the seed's lower 2 bytes into RNGEN's seed storage. 

POKE 8H162,Al!: POKE &H161,A2? 
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4270 GOTO 3500 'return to the menu 

4271 ' 

6500 !z============= CHANGE MODEL PARAMETERS SBRT / MPM zzzzzzzzzzzzzzzzzzzzzz 
4505 SCREEN O,1: COLOR 15,5,0: CLS 

4510 MSG1$ = “Enter a positive, nonzero number only." 

4515 MSG2$ = “Enter a positive value between 0 and 1." 

4520 LOCATE 1,6: PRINT “CHANGE MODEL PARAMETERS": PRINT BORDERS 


4525 LOCATE 4,1: INPUT “Enter Arrival Rate : ", LAMBDA! 

4530 IF LAMBDA'!«zO THEN LOCATE 23,1:PRINT MSG1$;:LOCATE 4,1:PRINT STRINGS$(40,32) 
:GOTO 4525 

4552 LOCATE 23,1: PRINT STRING$(40,232); 

4525 LOCATE 6,1: INPUT "Enter Service Rate for Server 1: ", Ul'! 

6560 IF Ul! <=0 THEN LOCATE 23,1: PRINT MSG1$3: LOCATE 6,1: PRINT STRINGS(40,32): 
GOTO 4535 


4542 LOCATE 23,1: PRINT STRING$(40,232); 

4565 LOCATE 8,1: INPUT “Enter Service Rate for Server 2: ", U2? 

4550 IF U2! «=O THEN LOCATE 23,1: PRINT MSG1$3: LOCATE 8,1: PRINT STRING$(40,32): 
GOTO 4545 

6555 LOCATE 23,1: PRINT STRING$(40,32); 

4560 LOCATE 10,1: INPUT "Enter Prob of Exit after Server 1: ',P1'! 

G565 IF ((P1! < 0) OR (P1?! » 1)) THEN LOCATE 23,1: PRINT MSG2$;: 
LOCATE 10,1: PRINT STRINGS$(40,32): GOTO 4560 

4567 LOCATE 23,1: PRINT STRING$(G40,32); 

4570 LOCATE 12,1: INPUT "Enter Prob of Exit after Server 2: ",P2! 

4575 IF ((P2* « 0) OR (P2! » 1)) THEN LOCATE 223,1: PRINT MSG2$;: 
LOCATE 12,1: PRINT STRINGS$(40,32): GOTO 4570 

4577 LOCATE 23,1: PRINT STRING$(40,32); 


4580 GOTO 120 'start simulation 
4581 ' 
GQ600 'zz2z2z2z2z2zzzz-z-- EXPLOSIVE QUEUE WARNING SBRT ¿Z MPM ss3zszszzzzzzz2z2zzzzzzz- 


G605 SCREEN 0,1: COLOR 14¢,3,0: WIDTH 40: CLS 

G610 LOCATE 6,5: PRINT “Your parameter selection results" 
G615 LOCATE 7,5: PRINT “in an Unstable (EXPLOSIVE) queue." 
4620 LOCATE 8,5: PRINT “The graphical presentation of the" 
G625 LOCATE 9,5: PRINT “network will still be accurate;" 
4630 LOCATE 10,5:PRINT "However the network statistics" 
4635 LOCATE 11,5:PRINT "will be INVALID, and will NOT be" 
G6G0 LOCATE 12,5:PRINT "displayed. See the user's manual" 
4661 LOCATE 13,5:PRINT "for an explanation." 

G645 LOCATE 15,10:PRINT “<C>ontinue with Simulation" 

4650 LOCATE 16,10:PRINT "<R>eturn to the Program MENU" 
G655 LOCATE 18,5: PRINT "ENTER your Selection" 

4660 GOSUB 29000 'Capture Keyboard input 

4665 AZ = INSTRI"CRer", RESPONSES ) 

4670 IF AZ = O THEN 4605 ELSE ON AZ GOTO 3720, 100, 3720, 100 
G671 ' 

G700 ‘============== DRAW & SAVE CIRCLES/ARROWS SBRT / MPM zzzzzzzzzzzzzzzzzz 
4701 ‘draw circles 

4702 CIRCLE (11,35),7,1:PAINT(11,35),1,1 

G706 Y=G : X=27 

4706 FOR I=Il TO 9 


4708 CIRCLE (X,35),7,2: PAINT(X,35),2,2 
4710 LOCATE 5,Y : PRINT RIGHTS(STRS(I),1) 
4712 X-X*16 

4714 Y=Y+2 

4716 NEXT 


4717 'draw arrows colored and umcolored 

G718 DRAW “c33 bm 05473 m4,4¢73h43 f43g43e4¢3 m135,473 hos £43 g4" 
4720 DRAW “c23 bm 0,633 m4,633h43443943e43 m1l35,63; há; f4; gá" 
4722 DRAW "c3, bm 159,67; m 271,473 há; f4; g4" 

4724 DRAW "c23 bm 159,633 m 271,63; há; f4; g4" 

4726 DRAW "c33 bm 303,793 m 15,793 £43 hGs eG" 

4728 DRAW "c3; bm 207,79; f4 h4; e4; bm 111,795; fé; hé; eé" 
4730 DRAW "c2; bm 303,95; m 15,95; f4; há; eG" 

4732 DRAW "c2; bm207,95; fé; hé; e4; bm 111,95; fé; hé; eé" 
4736 DRAM "c3; bm 25,47; m319,47; hé; fé; ga" 

4736 DRAW "c2; bm 295,633 m 319,635; hé; f4; g4" 

4738 DRAW “c3; bm 223,39; m 223,11; {άν Πάν σά" 

4740 DRAW "c2; bm 239,39; m 239,11; f4; há; g4" 
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4741 'save circles for future plotting 
4742 GET (4,28) - (18,42), CIRCO 
4744 GET (20,28) - (34,42), CIRC1 
4746 GET (36,28) - (50,42), CIRC?2 
4748 GET (52,28) - (66,42), CIRC3 
4750 GET (68,28) - (82,42), CIRC4 
4752 GET (84,28) - (98,42), CIRC5 
4754 GET (100,28) - (1145,42), CIRC6 
4756 GET (116,28) - (130,42), CIRC7 
4758 GET (132,28) - (146,42), CIRC8 
4760 GET (148,28) - (162,42), CIRC9 
4761 'save arrows for future plotting 
4762 GET(0,47) - (135,46), ARROMHI 
4764 GET(0,63) - (135,62), ARROHIC 
4766 GET(223,11) - (224,39), ARROWZ 
4768 GET(239,11) - (240,39), ARROW2C 
4770 GET(159,47) - (271,48), ARROWS 
4772 GET(159,63) - (271,64), ARROWSC 
4774 GET(2955,47) - (319,48), ARROWG 
4776 GET(295,63) - (319,64), ARRONGC 
4778 GET(15,79) - (303,80), ARROWS 
4780 GET(15,95) - (303,96), ARROWS5C 


4782 CLS 

4784 RETURN ‘return to line 3915 

4785 ' 

GB800 'zzzzzzzzzzzzzz-zzz MAIN SCREEN DRAHING SBRT / MPM 33a3z2z232z22z2z2z2z2z2zz2z2z22z2zz 


4801 'print network picture 

4802 LINE (0,8)-(319,8),2 : LINE (0, 10)-(319,10),2 
4804 LOCATE 1,14: PRINT "JACKSON QUEUE", 

4806 DRAW "c35; bm 0,43; m 135,435; Πάν f4; g4" 

4808 LINE (135,31)-(159,55),3,B 

4810 DRAW "c33 bm 159,435 m 271,433 hos £43 ge" 
4812 DRAW "c33 bm 167,433 m 167,15; f4; h4; g4" 
4814 LINE (271,31) - (295,55),3,B 

4816 DRAW "bm 2955433 m 319,43; hé; fé; gé" 

4818 DRAW "bm 3035433 m 303,763 h4s £43 e453 g^" 
4820 DRAWN "m 15,76; f4; h4; e4; g4" 

4822 DRAW "m 15,43; f4; h4; g4" 

4824 DRAN "bmé4,43; hé; fé; ga; bm 207,763 f4; hó; e4; bm 111,76; f4; hó; e4" 
4826 LOCATE 3,1: PRINT "lambda=" 

4828 LOCATE 3,23: PRINT "P1=" 

4830 LOCATE 3,34: PRINT "P2=" 

4832 LOCATE 7,7: PRINT "Ql=" 

4834 LOCATE 9,16: PRINT "ulz" 

4836 LOCATE 7,24: PRINT "q2z" 

4838 LOCATE 9,30: PRINT "u2=" 

4840 LOCATE 11,16: PRINT "TIME:" 

4841 ‘print state-verses time graphs 

4862 LINE (0,90) - (319,90),2 : LINE (0,92) - (319,92),2 
4844 LINE (157,92)-(157,199),2 

4846 LINE (23,185)-(23,104) 

4848 LINE (183,185)-(183,104) 

4850 FOR Y = 107 TO 179 STEP 8 


4852 LINE (19, Y) - (23,Y) 
4854 LINE (179,Y)-(183,Y) 
4856 NEXT 


4858 LINE (23,185)-(151,185) 
4860 LINE (183,185)-(311,185) 
4862 FOR X = 27 TO 151 STEP 4 


4864 IF (X=63) OR (X=103) OR (X=143) THEN GOTO 4870 
4866 LINE (X,185) - (X,188) 

4868 LINE (X+160, 185)-(X+160,188): GOTO 4874 

4870 LINE (X,185)-(X,190) 

4872 LINE (X*160,185)-(X*160,190) 

4874 NEXT 

4876 FOR COL = 1 TO 21 STEP 20: K = 10 

4878 FOR ROW = 15 TO 23 STEP 2 

4880 LOCATE ROW, COL 

4882 K=K-2 
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4884 
4886 
4888 
4890 
4892 
4896 
4896 
4898 
4900 
4902 
4904 
4906 
4908 
4910 
4912 
4916 
4916 
4918 
4920 
4922 
4926 
4926 
4928 
4930 
4932 
4934 
4936 
4937 
5000 
5010 
5030 
5031 
5032 
5033 
5036 
5050 
5080 
5081 
5085 
5088 
5090 
5095 
5100 


5110 
5115 
5120 
5125 
5126 
5127 
5130 
5135 
5140 
5144 
5145 
5146 
5147 
5148 
5149 
5150 
5151 
5152 
5155 
5156 
5160 
5161 
5162 
6000 
6010 


PRINT STRS(K) 


NEXT 


NEXT 


LOCATE 13,1: PRINT 
LOCATE 13,21: PRINT 
LOCATE 25,8: PRINT 
LOCATE 25,27: PRINT 


GET (158,93 )-(319,199), SERV2 
LINE (158,93 )-(319,199),0,BF 

= 0 THEN GOTO 4910 
24: PRINT "Unstable Queue" 


IF UNSTABLE 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
Y 
RETURN 


18, 
19,23: 
20,25: 
13,21: 
14,21: 
15,21: 
16,21: 
17,21: 
18,21: 
19,21: 
20,21: 
21,21: 
22,21: 
235,21: 
24,21: 


32: GOSUB 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
Ξ0000 


'return to 


"N( t)" 


"NÍ( t pM 


“SERVER 1°'5 


"SERVER 2"; 

'save server 2 state graph 
'erase server 2 state graph 
'stable do not print warning 


"Model Statistics" 
“Are Invalid": GOTO 4936 
“VAR ESTIMATE LIMIT" 
"R1 0.000" 
"R2 0.000" 
"SI 0.000" 
ES 0.000" 
"Ll 0.000" 
“LZ 0.000" 
v 0.000" 
"Wl 0.000" 
"W2 0.000" 
"W 0.000" 
uUZ 0.000"; 
‘print stat limits to screen 
line 3920 


PSHE EE EE EE EE EE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EHH Ete tet ee ett 


ΕΕ 1111111 


CLOCK MODULE Z CM *tttttttttttttttttttrtttàte 


TH HE EEE EEE EEE HES EE EE EE EEE EEE EEE EE EEE EEE EEE EE EEE EEE EEE EES EEE EE Hee t tte tt 


DEF SEG 


= 
-— 


NAT! 


= 
— 


GOSUB 25000 
GOSUB 26000 
LOCATE 7,27: PRINT USING FMT5$3 Q(2) 
LOCATE 7,10: PRINT USING FMT5$3 Q(1) 


IF ((Q(1) >= 250) OR (Q(2) >= 
IF UNSTABLE 1 THEN GOTO 5090 
IF (W>=N AND MODE 
FNMHIN!(NAT!, NSCT!(1), 
0 THEN GOSUB 7500 

= NAT! THEN NEXTEVNT = 1 ELSE IF CTIME! 


2 ELSE NEXTEVNT 


CTIME! 


IF CTIME! 
IF CTIME! 
NEXTEVNT 


NEXTEVNTTIME 


DEF SEG: POKE &H6A,0: 


&H1A00: CALL RNGEN(U' ) 
FNEXPONT*'(U*', LAMBDA'!) 
‘initialize exit(1) 
‘initialize exit(2) 


FNROUNDUPZ ( CTIME ! ) 


'generate a uniform(0,1) RV 
'initialize first exogenous arrival 
- exit after server 17: l=yes O=no 
- exit after server 27: l=yes O=no 


250)) THEN GOSUB 9600 ‘program crash 

‘skip updating statistics 
2) THEN W=O: GOSUB 9000 ‘update and print statistics 
NSCT!(2)) ‘select next event time 
'tie-breaking sbrt 
= NSCT'!(1) THEN 
'determine what next event is 
'round up to next integer 
'clear keyboard buffer 


3 


WHILE TIME < NEXTEVNTTIME 


TIME 


IF FAST 


TIME + 1 
1 THEN GOTO 5145 


'skip updating state graphs 


‘if time <= 32 skip screen shift sbrt, goto state diagram update sbrt 


IF TIME <= 
GOSUB 7000 


XX 


32 THEN GOSUB 6000: XX-2XX*4: FOR I 
‘screen shift subroutine if time >= 32 
147: GOSUB 6000 


1 TO 600: NEXT: GOTO 5145 


'state diagram update sbrt 


LOCATE 11,21: PRINT USING FMT1$3 TIME 


‘check keyboard 
INKEY$: 


AS 
IF 
IF 
IF 
IF 
IF 


(Αν 
( AŠ 
(A$ 
(A$ 
(AŠ 


DEF SEG: 


HEND 
W 


men 
a 
"d" 
“8 
ape 


"n # H H η 


buf 
ΑΦ 
AS 
AŠ 
AS 
EI 
ΑΦ 


1Ε 
OR 
OR 
OR 
OR 
OR 


POKE &H6A,0: 


SERVER 1 STATE DIAGRAM UPDATE SBRT / CM 
‘determines color of tic marks: 


fer for user input and act accordingly 
"" THEN GOTO 5152 ‘if no Keyboard input skip 
"S") THEN GOTO 100 
"P") THEN GOSUB 8000 : GOTO 5152 
= "D") THEN GOSUB 8500 : GOTO 5152 
"b") THEN GOSUB 8600 : GOTO 5152 
"f") THEN GOSUB 9500 
'clear keyboard buffer 


Hu n" wu 


W+l ‘event counter - determines when stats updated in line 5088 
ON NEXTEVNT GOTO 10000, 15000, 


‘arrival module, server 1 module or 
server 2 module 


20000 


l=blue 2=purple 
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6020 
6030 
6040 
6045 
6050 
6051 
6500 
6510 
6520 
6530 
6540 
6550 
6551 
7000 
7010 
7020 
7025 
7030 
7040 
7050 
7051 
7500 
7540 
7550 
7560 
7570 
7580 
1585 
7590 
7595 
1597 
7600 
7605 
7610 
7615 
7620 
8000 
8005 
8010 
8020 
8021 
8500 
8505 
8510 
8520 
8550 
8531 
8600 
8601 
8605 
8610 
8615 
8616 
9000 
9005 
9010 
9015 
9020 
9021 
9022 
9023 
9024 
9025 
9030 
9035 
9040 
9045 
9050 
9055 


NT1 = STAT(1) * Q(1) ‘total customers waiting and being served ἃὢ 1 
IF NTl > 9 THEN NTl1 = 9: Cl = 1 'can not plot more than 9 circles 
LINE (XX, 179 - (8*NT1)) - STEP(4,0), Cl ‘plot the tic mark 

IF MODE=1 THEN GOSUB 6500 ‘if server 2 state graph is displayed update 
RETURN ‘return to line 5130 or 5145 


.ΞΞΞΞΞΞΞΞΞ SERVER 2 STATE DIAGRAM UPDATE SBRT / CM ss=s=s================ 


NT2 = STAT(2) + Q(2) ‘total customers waiting and being served 232 
IF NT2 > 9 THEN NT2 = 9: C2 = 1 'can not plot more than 9 circles 
LINE (XX + 160, 179 - (8*NT2)) - STEP(G,0), C2 ‘plot the tic mark 
RETURN ‘return to line 5130 or 5145 


IS=5=7=s=2==z===22z222222 SCREEN SHIFT SBRI ο Giles === Sea eee ee 
GET (28, 103) - (155,183), SHIFT ‘capture server 1 state graph 


PUT (24,103), SHIFT, PSET ‘shift and re-plot the graph 

IF MODE = 2 THEN RETURN ‘if not displaying server 2 return 

GET (188,103) - (315,183), SHIFT ‘capture server 2 state graph 

PUT (184,103), SHIFT, PSET ‘shift and replot the graph 

RETURN ‘return to line 5140 

(================= TIE BREAKING SBRT / CM ============================= 
DEF SEG = &H1A00: CALL RNGEN(U!) 


IF NSCT!(1) = NSCT!(2) AND NSCT!(1) = NAT! GOTO 7590 ‘3-way tie 
IF NAT! = NSCT!(2) GOTO 7600 ‘2-way tie 
IF NAT! = NSCT!(1) GOTO 7610 ‘2-way tie 


IF U! <= .5 THEN CTIME! = NSCT!(1): NEXTEVNT = 2: GOTO 5110 ‘return clock 
CTIME! = NSCT!(2): NEXTEVNT = 3: GOTO 5110 
IF U! <= .333 THEN CTIME! = NAT!: NEXTEVNT = 1: GOTO 5110 ‘return to clock 


IF U? >= .333 AND U! <= .666 THEN CTIME! = NSCT!(1):NEXTEVNT=2: GOTO 5110 
CTIME! = NSCT!(2): NEXTEVNT = 3: GOTO 5110 ‘return to clock module 

IF Uf <= .5 THEN CTIME! = NAT!: NEXTEVNT = 1: GOTO 5110 ‘return to clock 
CTIME! = NSCT!(2): NEXTEVNT = 3: GOTO 5110 ‘return to clock module 

IF Ut <= .5 THEN CTIME! = NAT!: NEXTEVNT = 1: GOTO 5110 ‘return to clock 
CTIME! = NSCT!(1): NEXTEVNT = 2: GOTO 5110 ‘return to clock module 
!zzzzzzzz2zzzzzzzzzzzzz PAUSE SBRT / CM sszsz2zz2z2zzz2z2z22z22z2z222z2z2z2z2z2z2z2z2z2z2z2z2z2z 
'used in flashing of arrows 

DEF SEG: POKE &H6A,0: 


AS = INKEYS: IF AS="c" OR AS$="C" THEN RETURN ELSE GOTO 8020 


'z========= CHANGE DISPLAY CONFIGURATION SBRI PZ CH τ τ τι; 
'replaces server 2 state diagram with stats table or vice/versa 

IF MODE = 1 THEN PUT(158,93), STATISTICS, PSET: MODE = 2: RETURN 

PUT (158,93), SERV2, PSET: MODE = 1 

RETURN ‘return to line 5149 


tz================ BEEP SUPPRESSION SBRT / CM szsszszszzszzszzzzzzzzzzzzzzzzz 
‘turns the audible beeps on and off 
IF DUR1 = 3 THEN DUR1 = O ELSE DURI 
IF DUR2 = 2 THEN DUR2 = O ELSE DUR2 
RETURN ‘return to line 5150 


3 
2 


*‘sss=ssss= UPDATE STATISTICAL ESTIMATES SBRT / CM s====================s== 
TA1 = NA1 + FB1 
R1! = TA1 Z CTIME! 


TA2 Z CTIME! 
IF CUMSERV1! = O THEN SR1! = 0: SR2! = 0: GOTO 902% 
SR1! = (El + TA2)/CUMSERV1! ‘estimate of service rate for server 1 
IF CUMSERV2! = 0 THEN SR2 Ξ 0: GOTO 9024 
SR2! = (E2 + FB1)J/CUMSERV2! ‘estimate of service rate for server 2 
R! = NAIZCTIME! 'estimate of exogenous arrival rate 
S1!=O: S2!=O: ET= El + E2: IF ET = O THEN GOTO 9056 
FOR I = 1 TO 10 
S1! S1! + IwSERVBY1(I) 
S2! S2' + IXSERVBY2(I) 
NEXT 
S1! - S1'/ET 
S2! = S2'/ET 
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9056 
9060 
9062 
9065 
9070 
9075 
9077 
9080 
9085 
9086 
9087 
9088 
9090 
9200 
9210 
9211 
9500 
9501 
9505 
9510 
3515 
9520 
9525 
9530 
9535 
9540 
9545 
9600 
9605 
9610 
9615 
9620 
9621 
9625 
9626 


10000 
10010 
10020 
10030 
10035 
10037 
10040 
10042 
10045 
10050 
10060 
10063 
10065 
10070 
10080 
10090 
10091 
11000 
11005 
11010 


11020 
11030 
11021 
12000 
12010 
12020 
12030 


IF (SRI! «= RI!) THEN LI! = 0: GOTO 9062 

LI! = RI!Z(SRI! - RI!) 

IF (SR2! «= R2!) THEN L2! = 0: GOTO 9070 

L2! = RZ2IZ(SRZ2! - R2!) 

L! = LI! + L2! 

W1! = L1!/R1! 

IF (R2! = 0Ο) THEN W2!=0: GOTO 9085 

W2! - L2'/R2'! 

W! = LIZR! 

IF SRI! = O THEN Z! = 0: GOTO 9200 

IF (P1! = 1) THEN Z! = (RI'!/SRI!): GOTO 9200 
IF SR2! = O THEN Z! = 0: GOTO 9200 

Z! = ABS((R1!⁄7SR1!)X(1 - (R2!/SR2!))) 

Y=23: GOSUB 30000 ‘plot estimates to screen 
RETURN ‘return to line 5090 


!zzzzzzzzzzzzz- GRAPHICS SUPPRESSION SBRT / CM zzszzzzzzzzzzzzzzzzzzzzzzzz 
‘1f fast mode is selected no graphics are displayed 
IF FAST = O THEN FAST = 1 ELSE FAST = 0: N=1 
IF FAST = 1 THEN GOTO 9515 ELSE RETURN 
N = 10 ‘when fast mode selected, stats updated after every N events 
LINE (140,36) - (154,50),0,BF ‘erase customers and state tic marks 
LINE (4,28) - (130,42),0,BF 
LINE (276,36) - (290,50),0,BF 
LINE (172,28) - (266,42),0,BF 
LINE (24,105) - (151,183),0,BF 
RETURN ‘return to line 5152 
!"Zzzzzzzzzzzzzzzz PROGRAM CRASH SBRT / CM zzzzszzzzzzzzzzzzzzzzzzzzzzzzzz 
SCREEN 0,1: COLOR 14,3,0: WIDTH 40: CLS 
LOCATE 7,3: PRINT "Program execution HALTED. Queue" 
LOCATE 8,3: PRINT “length exceeds program capabilities." 
LOCATE 9,3: PRINT “Press any Key to return to Main" 
LOCATE 10,3: PRINT "Program menu."; 
GOSUB 29000: GOTO 100 'display program menu 
t6 6G GAB B  . B B. t.  . | t. + tt tr tt tit tt tr tt tt tt tt tt tr tt tt tt+r tt tr th tt Ht tt +K tt te++ 
!++++++++++++++++++ ARRIVAL MODULE / AM +Y+++++++++++++++++++++++++++ 
T+T+TT+Ti+TÑiliitiiiliT+iii+iiiii+i+i+iTTTTTTTTTTTTTTTTTTTTTTTTTTT71TV74+ +++++++4+++++++++ 
STRG1$ = "0" + STRG1$ 'O corresponds to new arrival, add to queue strg 
ΝΑΙ = NA1 + 1 
IF FAST = 1 THEN GOTO 10042 'skip plotting routine 
IF Q(1)<8 THEN GOSUB 11000 ELSE GOSUB 12000: PUT(0,42),ARROW1,PSET 
DEF SEG = &H1A00: CALL RNGEN(UE ) 
AINC! = FNEXPONT!(U!, LAMBDA!) ‘compute next inter-arrival time 
NAT! = AINC! + CTIME! 'compute new exogenous arrival time 
IF STAT(1) = 1 THEN Q(1) = Q(1) + 1: GOTO 5080 'return to clock module 
CALL RNGEN(U! ) 
S1INC! = FNEXPONT!(U!? ,U1? ) ‘compute new server 1 service duration 
NSCT!(I) = SIINC! + CTIME! ‘compute next service completion time for 1 
STAT(1) 1 ‘change status of server 1 to busy 
GOTO 5081 ‘return to the clock module 
' 


A ooo ARRIVAL PLOTTING SBRT / AM zzzzzzzzzzzzzzzzzzzzzzzzzz 

GOSUB 12000 ‘flash arrow 1 

IF STAT(1) = O THEN PUT(1460,36), CIRCO:GOSUB 28000: PUT(0,42), ARROWI,PSET: 
RETURN ‘return to line 10040 

PUT (116-Q(1)*16,28), CIRCO:GOSUB 28000: PUT (0,42), ARROW1, PSET 


RETURN ‘return to line 10040 

!zzzzzzzzzzzzzz2ez FLASH ARROW l1 SBRT / AM ============================s= 
PUT(0,42), ARROWIC, PSET 'flash purple arrow 

FOR I = 1 TO 500 

NEXT 


12035 SOUND 500,DURI 


12040 
12041 
15000 


RETURN ‘return to line 10040 or 11010 


696A B RR 9 9. 9. 9. 9. 9. 4. 9. 4 9. 9 94 94 9. 94 9. 9 9 9. 9. 9. 9. it 4. 9. 4 4 4 4 9. 9. 4. 4 49 4 4 9 9 9 9 9 9 tt tt tt tit tt tt tr tét tt tete N 


15010 '++++++++++++++++ SERVICE CENTER 1 HODULE / SCIM +++++++++++++++++++++ 
15020 '+++++++++++++++++++++++++++4i+0i -d++++++++++++++++RT+T+++T+++++++++++++++++++++ 
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15025 
15030 
15035 
15040 
15050 
15060 
15070 
15072 
15075 
15080 
15090 
15099 
15100 
15110 
15113 
15115 
15120 
15130 
15160 
15161 
17000 
17010 
17020 
17030 
17035 
17040 
17050 
17055 
17060 
17070 


17080 
17090 
17095 
17099 
17100 
17110 
17115 
17120 
17130 
17140 
17150 
17156 
17155 
17160 
17170 
17175 
17180 
17181 
18000 
18010 
18015 
18020 
18030 
18040 
18050 
18051 
18100 
18110 
18120 
18130 
18135 
18140 
18141 
18500 
18501 
18510 
18520 
18530 


CUMSERV1! = CUMSERV1! + SLINC! ‘update amount time server 2 busy 
TRANSFER1S$ = RIGHTS(STRG1$,1) ‘move customer from Ql to Q2 
IF EXIT(1) = O THEN STRG2$ = TRANSFERIS + STRG2$: TA2 = TA2+1 
STRG1$ - LEFT$( STRG1$, LEN(STRG1$) - 1) 
IF FAST = 0 THEN GOSUB 17000  'server 1 plotting sbrt 
IF EXIT(1) 1 THEN El = El + 1: GOSUB 18500: GOTO 15100 ‘estimate 51 
IF STAT(2) 1 THEN Q(2) = Q(2) + 1 : GOTO 15100 
DEF SEG &H1A00: CALL RNGEN(U!) 
SZINC! = FNEXPONT'!(U!, U2!) |'compute new server 2 service duration 
NSCT!(2) = S2INC! + CTIME! 'compute new service completion time 
STAT(2) = 1 "change status for server 2 to busy 
‘make server 1 completion impossible event 
IF Q(1) » O THEN NSCT!(1) = 1E+3]l: STAT(1) = 0: GOTO 15130 


Q(1) -Q(1) - 1 
DEF SEG = &H1A00: CALL RNGEN(U! ) 
S1INC! = FNEXPONT!(U!, U1!) ‘compute new server 1 service duration 
NSCT!(1) = SLINC! + CTIME! ‘compute new service completion time 
GOSUB 25000 ‘generate new exit(1) indicator variable 
GOTO 5050 ‘return to clock module 
Emm mm mmm mmm mmo ms rm SERVER 1 PLOTTING SBRT 7 SCIM ΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞ-Ξ 
GET( 140,36) - (154,50), TEMPCIRC ‘capture circle a server l 
LINE(140,36) - (154,50),0,BF ‘erase circle d server l 
IF EXIT(1) = 1 THEN GOSUB 18000: GOTO 17060 ‘flash arrow 2 


GOSUB 18100 ‘flash arrow 3 

ΙΕ STAT(2) = O THEN PUT(276,36), TEMPCIRC : GOTO 17055 ‘put at server 2 
IF Q(2)< 6 THEN PUT(252-(Q(2)*16),28), TEMPCIRC ‘put in server 2 queue 
GOSUB 28000 :PUT(159,43), ARROW3, PSET ‘plot white arrow 3 


IF Q(1) = O THEN RETURN ‘return to line 15060 
GET(116,28)-(130,42), TEMPCIRC : LINE (116,28)-(130,42),0,BF: 
PUT(140,36), TEMPCIRC ‘put next customer in server l box 
IF Q(1) = 1 THEN RETURN ‘return to line 15060 
x = 100 
IF Q(1)<= 8 THEN K = (Q(1) -1) ELSE K = 7 


'shift all remaining customers to the right one 
FOR J = 1 TO K 
GET(X,28) - (X+14,42), TEMPCIRC 
LINE (X,28) -(X*164,42),0,BF 
PUT(X+16,28), TEMPCIRC 
x = X-16 
NEXT 
IF Q(1)<9 THEN RETURN ‘return to line 15060 


'go to overflow queue and bring next customer into visible queue 
LOOPCOUNT1$ z MIDS(STRG1$,(LENCSTRG1$)-8),1) 

X=4 : Y=28 

INDEX * VALCLOOPCOUNT1$ ) +1 

GOSUB 27000  'bring circle in from queue overflow 

RETURN ‘return to lne 15060 


mE mmm FLASH ARROM 2 SBRT / SClM szzzzzzzzzzzzzzzzzzzzzzzz 
PUT (167,15), ARROW2C, PSET ‘flash purple arrow 2 


SOUND 3000,DUR2 

FOR I=1 TO 1000 

NEXT 

PUT(167,15), ARROW2, PSET ‘plot white arrow 2 

RETURN ‘return to line 17030 

ae ee ce eee ee FLASH ARROW 3 SBRT 7/7 SCIM zszzzszzzzzzzzzzzzzzzzzzzzz 
PUT (159,43), ARROW3C, PSET ‘flash purple arrow 3 


FOR I=1 TO 500 

NEXT 

SOUND 500,DUR1 

RETURN ‘return to line 17030 

"SSSS==5=== FREQUENCY OF SERVICE BY SERVER 1 SBRT / SCIM =============== 
"this sbrt used in calculation of S1 

FOR I = O TO 9 

IF VALCTRANSFER1$) = I THEN SERVBY1(I+1) = SERVBY1(I+1) + 1: GOTO 18540 
NEXT 
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18535 
18540 
18550 
18560 
18565 
18566 
20000 
20010 
20020 
20025 
20027 
20030 
20045 
20050 


20060 
20070 
20080 
20090 
20100 
20102 
20105 
20110 
2010 
20129 
20130 
20140 
20163 
20165 
20150 
20160 
20170 
20171 
22000 
22010 
22020 
22030 
22035 
22040 
22050 
22055 
22060 
22069 
22070 


22080 
22090 
22100 
22109 
22110 
22120 
22130 
22140 
22150 
22160 
22170 
22179 
22180 
22190 
22195 
22200 
22205 
22206 
23000 
23010 
23015 
23020 
23030 
23040 


IF VAL(TRANSFER1$) = 
FOR I = 1 TO 9 

IF VAL(TRANSFERI$) = 
NEXT 
RETURN 


O THEN RETURN 


I THEN SERVBY2(I) - SERVBY2(I) + 1: RETURN 


'return to line 15060 


E AB B B BG GAB GA VA. .,  ., 4.4. 4. t... -. 4€. 4 4.  . . 4... 4. 4. 46 6 6 B G A À Á AÀ VÀ V À BA ÀBA AG Ó GS 9 V V B t. t.t tB6A 

"4¢4¢444444¢444¢+4+ SERVICE CENTER 2 MODULE 7/7 SCZ2M ++++++++++++++++++++++++ 

EEE EEEE EEE EE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EEE EE 9 9 9 9 4 4 94 9. 444 

CUMSERV2Z! = CUMSERV2Z! + S2INC! ‘update amount of time server 2 busy 

TRANSFER2$ = RIGHTS(STRG2$S$,1) 

IF EXIT(2) = 1 THEN E2 = E2+1: GOSUB 23500: GOTO 20070 

FBI = FBI + 1 ‘increment feedback at server 1 counter 

IF (VAL(TRANSFER2$)«9) THEN TRANSFER2$S$ =RIGHTS( STRS( 
ELSE TRANSFERZ$ = "9e 

STRG1$ = TRANSFER2Z$ + STRG1$ 

STRG2$ - LEFTS(STRG2$, LEN(STRG2$) - 1) 

IF FAST = O THEN GOSUB 22000 ‘server 2 plotting subrt 

IF EXIT(2) = 1 THEN GOTO 20130 

IF STAT(1) = 1 THEN Q(1) = Q(1) 41 : 

DEF SEG = &H1A00: CALL RNGEN(U! ) 

SIINC! = FNEXPONT'!(U!, U1!) 'compute new server l service duration 

NSCT'(1) » SIINC! « CTIME'! 'compute new service completion time 

STAT(1) = 1 'change status of server 1 to busy 

‘make service completion at server 2 impossible event 

IF Q(2) = O THEN NSCT!#(2) = 1&4+31: STAT(2) = 0: GOTO 20160 

Q(2) = Q(2) - 1 

DEF SEG = &H1A00: CALL RNGEN(U!) 

SZINC! := FNEXPONT'(U!, U2!) 'compute new server 2 service duration 

NSCT!(2) = SZINC! + CTIME! 'compute new service completion time 

GOSUB 26000 'generate new exit(2) indicator variable 

GOTO 5050 'return to clock module 


'estimate S2 


VAL(TRANSFER2$)+1),1) 


GOTO 20130 


SERVER 2 PLOTTING SBRT / SC2M 

LINE (276,36) - (290,50),0,BF  'erase circle 9 server 2 

IF EXIT(2) * 1 THEN GOSUB 23000: GOTO 22060 ‘flash arrow 4 

INDEX = VAL(TRANSFERZ$) +1 ‘increment feedback counter 

GOSUB 23100 ‘flash arrow 5 

IF STAT(1) = O THEN Xz140 : Y-36 : GOSUB 27000 : GOTO 22055 

IF Q(1) < =7 THEN X = 116 - Q(1) * 16 : Y=28 : GOSUB 27000 

GOSUB 28000: PUT(15,76), ARROWS, PSET ‘plot white arrow 5 

IF Q(2) = 0 THEN RETURN ‘return to line 20090 

‘erase ball at server 1 and place at server 2 

GET( 252,28) - (266,42), TEMPCIRC LINE( 252,28 )-(266,42),0,BF: 
PUT(276,36), TEMPCIRC 


IF Qi2) = 1 THEN RETURN ‘return to line 20090 
%=236 
IF Q(2) <= 6 THEN K = (Q(2)-1) ELSE K = 5 


‘shift all remaining customers one to the right 
FOR J = 1 TO K 
GET(X,28) - (X+14¢,42), TEMPCIRC 
LINE(X,28) - (X*14,32), 0, BF 
PUT(X*16,28), TEMPCIRC 
X=X-16 
NEXT 
IF Q(2) <= 6 THEN RETURN 
‘go to overflow queue and bring next customer into visible queue 
LOOPCOUNT2$ = MIDS(STRG2$S$, (LEN(ISTRG2$)-6),1) 
x = 172 : Y = 28 
INDEX = VAL(LOOPCOUNT2$) + 1 
GOSUB 27000 ‘bring circle in from queue overflow 
RETURN ‘return to line 20090 
f 
t================ FLASH ARRON 4 SBRT / SC2M =s===s=s=s=======s=sx=s=sxz=s=====s===s=s=zs= 
PUT(295,43), ARROW4C, PSET ‘flash purple arrow 4 
SOUND 3000, DUR2 
FOR I = 1 TO 1000 
NEXT 
PUT (295,43), ARROWG, PSET ‘plot white arrow 4 
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23050 
23051 
23100 
23110 
23120 
23130 
25135 
23140 
23141 
23500 
23501 
23510 
23520 


23530 
24990 
24991 
24992 
24993 
24996 
25000 
25004 
25005 
25010 
25015 
25016 
26000 
26004 
26005 
26010 
26015 
26016 
27000 
27010 


27020 
27030 
27040 
27050 
27060 
27070 
27080 
27090 
27100 
27110 
27111 
28000 
28010 
28020 
28030 
28031 
29000 
29005 
29010 
29020 
29021 
30000 
30008 
30009 
30010 
30015 
30020 
30025 
30030 
30035 
30040 
30045 
30050 
30055 


RETURN ‘return to line 22020 
ΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞ FLASH ARROW 5 SBRT / SC2M 

PUT( 15,76), ARROWSC, PSET ‘flash purple arrow 5 
FOR I=1 TO 500 

NEXT 

SOUND 500,DUR1 

RETURN ‘return to line 22040 


FREQUENCY OF SERVICE BY SERVER 2 SBRT / SC2M 

'this sbrt used in calculation of S2 

FOR I=0 TO 9 

IF VAL(TRANSFER2$) * I THEN SERVBY2(I*1) - 
SERVBY1(I+1) = SERVBY1L(I+1) + 1: RETURN 

NEXT 


SERVBY2(I+l) + l: 


! 3€3€ € € € 9€3€ 3€ 9€ € 9€ X 9€ € 9€ 3€ € 9€ 9€ 9€ 3€ 9€ € 9€3€ 9€ € € 9€ C € 3€ 39€ 3€ 3€ € JE JE JE JE IE IE FE FETE JE JE JE JE JE JEJEJE JE JE JE JEJE JE IE IE JE IEIC II ICE 


' 333009. SUBROUTINES USED BY MORE THAN ONE MODULE 363636363636X3XX33633X 
! X E JEJE EJE JEJE E OC CE EXC € € JE JE JE JE € € JE JE JE JE JE PE JE PE JE JE JEJE JE HE JE JE JE JE JE JE JE JE JEFE III III E JE JE JE JE IEE EE 


PROB EXIT AFTER SERVER 1 SBRT / CM, 
‘exit after service at server 1: l=yes O=no 
DEF SEG = &H1A00: CALL RNGEN(U! ) 

IF Uf <= Pl! THEN EXIT(1) = 1 ELSE EXIT(1) = O 


SC1M, 


RETURN ‘return to line 5034 or 15140 

‘s=s=s===== PROB EXIT AFTER SERVER 2 SBRT / CM, SC2M zzzzzzzzzzzzzzzzzzzzz 
‘exit after service at server 2: l=yes O=no 

DEF SEG = &H1A00: CALL RNGEN(U!? ) 

IF Ufc= P2! THEN EXIT(2) = 1 ELSE EXIT(2) = 0 

RETURN ‘return to line 5050 or 20170 

!zzzzzzzzzzzzzzz QUEUE OVERFLOM SBRT / SC1M, πο) == eee ee 
ON INDEX GOTO 27020, 27030, 27040, 27050, 27060, 27070, 27080, 


27090, 27100, 27110 
PUT(X,Y), CIRCO : RETURN ‘return to line 17180 or 22040 or 22205 
PUT(X,Y), CIRC1 : RETURN 
PUT(X,Y), CIRC2 : RETURN 
PUT(X,Y), CIRC3 : RETURN 
PUT(X,Y), CIRC4 : RETURN 
PUT(X,Y), CIRC5 : RETURN 
PUT(X,Y), CIRC6 : RETURN 
PUT(X,Y), CIRC7 : RETURN 
PUT(X,Y), CIRC8 : RETURN 
PUT(X,Y), CIRC9 : RETURN 
!zZzzzzzzzzzzzzz-cz DELAY SBRT / SCIM, SC2M szszzzzzzz2z2zzzzzzzzzzzzzzzzzz 
FOR I = 1 TO 500 
NEXT 
RETURN 
ΓΞΞΞΞΞΞΞΞΞΞΞΞ CAPTURE KEYBOARD INPUT SBRT Z MPM, CM ΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞΞ 
DEF SEG: POKE 8H6A,0: 
RESPONSE$ = INKEY$: IF RESPONSES = "" THEN 29010 
RETURN 
!zzzzzzzzzzzzzzc PRINT STATISTICS SBRT 7 MPM, CM ======================= 
'sbrt prevents displaying number that exceed formats 
IF (.R1* « .001 AND R1?!»0) OR (R1?! » 9999! )) THEN R1! = 9999.999 
LOCATE 16,Y: PRINT USING FMT2$; RI! 
IF ((R2* « .001 AND R2'!»0) OR (R2! » 9999!)) THEN R2! - 9999.999 
LOCATE 15,Y: PRINT USING FMT2$; R2! 
IF ((SI! « .OOI AND SI!»O) OR (SI! » 9999!)) THEN SI! = 9999.999 
LOCATE 16,Y: PRINT USING FMT2$3 Sl! 
IF ((S2* « .001 AND S2!'>0) OR (S2! > 9999!)) THEN S2! = 9999.999 
LOCATE 17,Y: PRINT USING FMT2$; S2! 
IF ((L1! < .001 AND L1!>0) OR (LI! > 9999!)) THEN LI! = 9999.999 
LOCATE 18,Y: PRINT USING FMT2$; L1'! 
IF ((L2! < .001 AND L2!>0) OR (L2?! » 9999? )) THEN L2! = 9999.999 
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30060 
30065 
30070 
30075 
30080 
30085 
30090 
30095 
30100 
30105 
30110 
30130 


LOCATE 19,Y: PRINT 
IF ((L! < .001 AND 
LOCATE 20,Y: PRINT 


IF (CWI! « .001 AND W1!>0) OR (W1! 


LOCATE 21,Y: PRINT 


IF ((N2* « .001 AND W2!'>0) OR (W2! 


LOCATE 22,Y: PRINT 
IF ((W! < .001 AND 
LOCATE 23,Y: PRINT 
IF ((Z! « .001 AND 
LOCATE 26,Y: PRINT 
RETURN ‘return to 


USING FMT2$3 L2! 


L'>0) OR (LY > 99991)) THEN L! 


USING FMT2$3 L! 
USING FMT2$; WI! 


USING FMT2$3 W2! 


W!>0) OR CW! > 99991)) THEN H! 


USING FMT2$; W!; 


Z1>0) OR (Z2! » 9999?)) THEN Z! 


USING FMT2$3 Z'!; 
line 4936 or 9210 


> 9999! )) THEN 


» 9999!)) THEN 
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- 9999,999 


W1! 2999. 999 


We 9999; 999 


9999.999 


9999.999 
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